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Abstract 

We argue that the natural functions for describing the multi-Regge limit of six- 



\Q ! gluon scattering in planar = 4 super Yang-Mills theory are the single-valued 

harmonic polylogarithmic functions introduced by Brown. These functions depend 
on a single complex variable and its conjugate, (w, w*). Using these functions, and 
■ formulas due to Fadin, Lipatov and Prygarin, we determine the six-gluon MHV 

remainder function in the leading-logarithmic approximation (LLA) in this limit 
through ten loops, and the next-to-LLA (NLLA) terms through nine loops. In sep- 
arate work, we have determined the symbol of the four-loop remainder function for 
general kinematics, up to 113 constants. Taking its multi-Regge limit and match- 
ing to our four-loop LLA and NLLA results, we fix all but one of the constants 
that survive in this limit. The multi-Regge limit factorizes in the variables (u, n) 
which are related to (w,w*) by a Fourier- Mellin transform. We can transform 
the single-valued harmonic polylogarithms to functions of (u, n) that incorporate 
harmonic sums, systematically through transcendental weight six. Combining this 
information with the four-loop results, we determine the eigenvalues of the BFKL 
kernel in the adjoint representation to NNLLA accuracy, and the MHV product 
of impact factors to N 3 LLA accuracy, up to constants representing beyond-the- 
symbol terms and the one symbol-level constant. Remarkably, only derivatives of 
the polygamma function enter these results. Finally, the LLA approximation to 
the six-gluon NMHV amplitude is evaluated through ten loops. 



1 Introduction 



Enormous progress has taken place recently in unraveling the properties of relativistic scattering 
amplitudes in four-dimensional gauge theories and gravity. Perhaps the most intriguing devel- 
opments have been in maximally supersymmetric M = 4 Yang-Mills theory, in the planar limit 
of a large number of colors. Many lines of evidence suggest that it should be possible to solve 
for the scattering amplitudes in this theory to all orders in perturbation theory. There are also 
semi-classical results based on the AdS/CFT duality to match to at strong coupling [I]. The 
scattering amplitudes in the planar theory can be expressed in terms of a set of dual (or region) 
variables , which are related to the usual external momentum four- vectors fcf by fcj = Xi —Xi+i. 
Remarkably, the planar M = 4 super- Yang-Mills amplitudes are governed by a dual conformal 
symmetry acting on the Xi [1] El El El El El E] • This symmetry can be extended to a dual super- 
conformal symmetry [8], which acts on supermultiplets of amplitudes that are packaged together 
by using an N = 4 on-shell superfield and associated Grassmann coordinates [H [TO], [TTJ [12] . 

Due to infrared divergences, amplitudes are not invariant under dual conformal transforma- 
tions. Rather, there is an anomaly, which was first understood in terms of polygonal Wilson 
loops rather than amplitudes [7]. (For such Wilson loops the anomaly is ultraviolet in nature.) 
A solution to the anomalous Ward identity for maximally-helicity violating (MHV) amplitudes 
is to write them in terms of the BDS ansatz [TBI. 



^MHV = A BBS x exp(/U (L1) 



where R n is the so-called remainder function [T4l[T5], which is fully dual-conformally invariant. 

For the four- and five-gluon scattering amplitudes, the only dual-conformally invariant func- 
tions are constants, and because of this fact the BDS ansatz is exact and the remainder function 
vanishes to all loop orders, = R$ = 0. For six-gluon amplitudes, dual conformal invariance 
restricts the functional dependence to have the form i? 6 (wi, w 2 , W3), where the Ui are the unique 
invariant cross ratios constructed from distances x\a in the dual space: 

22 22 22 

X 13 X 46 _ S12S45 _ X 24 ^X 15 _ S23S56 _ £35^26 _ 'S34S6I . > 

U\ — — 2~ — ) u 2 — — 2~ ~ ' U 3 ~ ~~2 2~~ ~~ • l 1 -^/ 

x 14 x 36 S 123S345 x 25 x 14 S 234S456 X 36 X 25 S 345S561 

The need for a nonzero remainder function R n for Wilson loops was first indicated by the strong- 
coupling behavior of polygonal loops corresponding to amplitudes with a large number of gluons 
n [6]. At the six-point level, investigation of the multi-Regge limits of 2 — > 4 gluon scattering 
amplitudes led to the conclusion that Rq must be nonvanishing at two loops [16]. Numerical 

(2) 

evidence was found soon thereafter for a nonvanishing two-loop coefficient R 6 for generic non- 
singular kinematics [13], in agreement with the numerical values found simultaneously for the 
corresponding hexagonal Wilson loop |15j . 

Based on the Wilson line representation JT5] . and using dual conformal invariance to take a 
quasi-multi-Regge limit and simplify the integrals, an analytic result for Rq was derived [H] HE] 
in terms of Goncharov's multiple polylogarithms [19] . Making use of properties of the symbol [20], 
|2T| [22] [23] 124"] associated with iterated integrals, the analytic result for R 6 was then simplified 
to just a few lines of classical polylogarithms [23] . 
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A powerful constraint on the structure of the remainder function at higher loop order is pro- 
vided by the operator product expansion (OPE) for polygonal Wilson loops [251 [26J EZ]- At 
three loops, this constraint, together with symmetries, collinear vanishing, and an assumption 

(3) 

about the final entry of the symbol, can be used to determine the symbol of R 6 up to just two 
constant parameters [28]. Another powerful technique for determining the remainder function is 
to exploit an infinite-dimensional Yangian invariance [291 US] which includes the dual supercon- 
formal generators. These symmetries are anomalous at the loop level (or alternatively one can 
say that the algebra has to be deformed) [3T] ■ However, the symmetries imply a first order linear 
differential equation for the £-loop n-point amplitude, and the anomaly dictates the inhomoge- 
nous term in the differential equation, in terms of an integral over an [I — l)-loop (n + l)-point 
amplitude [32], [33] . Using this differential equation, a number of interesting results were obtained 
in ref. [33]. In particular, the result for the symbol of Rq found in ref. [28] was recovered and 
the two previously-undetermined constants were fixed. 

In principle, the method of refs. [521 E2] works to arbitrary loop order. However, it requires 
knowing lower-loop amplitudes with an increasing number of external legs, for which the number 
of kinematical variables (the dual conformal cross ratios) steadily increases. Although the symbol 
of the two-loop remainder function R n is known for arbitrary n [31], the same is not true of 
the three-loop seven-point remainder function, which would feed into the four-loop six-point 
remainder function — one of the subjects of this paper. 

In this article, we focus on features of the six-point kinematics that allow us to push directly 
to higher loop orders for this amplitude, without having to solve for amplitudes with more legs. 
In fact, most of our paper is concerned with a special limit of the kinematics in which we can 
make even more progress: multi-Regge kinematics (MRK), a limit which has already received 
considerable attention in the context of Af = 4 super- Yang-Mills theory [161 [23 ESI [361 E3 EU ESI 
EH SI1 H21 H3J • In the MRK limit of 2 — > 4 gluon scattering, the four outgoing gluons are widely- 
spaced in rapidity. In other words, two of the four gluons are emitted far forward, with almost 
the same energies and directions of the two incoming gluons. The other two outgoing gluons are 
also well-separated from each other, and have smaller energies than the two far-forward gluons. 

The MHV amplitude possesses a unique limit of this type. For definiteness, we will take 
legs 3 and 6 to be incoming, legs 1 and 2 to be the far-forward outgoing gluons, and legs 4 and 
5 to be the other two outgoing gluons. Neglecting power-suppressed terms, helicity must be 
conserved along the high-energy lines. In the usual all-outgoing convention for labeling helicities, 

the helicity configuration can be taken to be (+- 1 H — )• For generic 2 — > 4 scattering in 

four dimensions there are eight kinematic variables. Dual conformal invariance reduces the eight 
variables down to just the three dual conformal cross ratios Ui. Taking the multi-Regge limit 
essentially reduces the amplitude to a function of just two variables, w and w*, which turn out 
to be the complex conjugates of each other. 

We will argue that the function space relevant for this limit has been completely characterized 
by Brown [H]. We call the functions single-valued harmonic polylogarithms (SVHPLs). They 
are built from the analytic functions of a single complex variable that are known as harmonic 
polylogarithms (HPLs) in the physics literature [IS]. These functions have branch cuts at w = 
and w — — 1. However, bilinear combinations of HPLs in w and in w* can be constructed [H] 
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to cancel the branch cuts, so that the resulting functions are single-valued in the (w,w*) plane. 
The single- valued property matches perfectly a physical constraint on the remainder function in 
the multi-Regge limit. SVHPLs, like HPLs, are equipped with an integer transcendental weight. 
The required weight increases with the loop order. However, at any given weight there is only a 
finite-dimensional vector space of available functions. Thus, once we have identified the proper 
function space, the problem of solving for the remainder function in MRK reduces simply to 
determining a set of rational numbers, namely the coefficients multiplying the allowed SVHPLs 
at a given weight. 

In order to further appreciate the simplicity of the multi-Regge limit, we recall that for generic 
six-point kinematics there are nine possible choices for the entries in the symbol for the remainder 
function Re(ui, u 2 , u 3 ) [23| 1251: 



{ill, «2, «3, 1 - «1, 1 - U2, 1 - "3, Vl, 2/2, 1/3} , (1.3) 

where 

Vi = (1.4) 
m - z_ 

-l + ui + u 2 + u 3 ± A 

*± = 2 ' ( ^ 

A = (1 - m - u 2 - w 3 ) 2 - 4uiu 2 u 3 . (1.6) 

The first entry of the symbol is actually restricted to the set {111,112,113} due to the location of 
the amplitude's branch cuts [27]; the integrability of the symbol restricts the second entry to the 
set {ui, 1 —U{} [271128]: and a "final-entry condition" [28| IM] implies that there are only six, not 
nine, possibilities for the last entry. However, the remaining entries are unrestricted. The large 
number of possible entries, and the fact that the yi variables are defined in terms of square-root 
functions of the cross ratios (although the Ui can be written as rational functions of the [28J), 
complicates the task of identifying the proper function space for this problem. 

So in this paper we will solve a simpler problem. The MRK limit consists of taking one of the 
Mi, say Ui, to unity, and letting the other two cross ratios vanish at the same rate that U\ — > 1: 
u 2 ~ x(l — ui) and u 3 m y(l—Ui) for two fixed variables x and y. To reach the Minkowski version 
of the MRK limit, which is relevant for 2 — > 4 scattering, it is necessary to analytically continue 
Mi from the Euclidean region according to u\ — > e~ 2m \ui\, before taking this limit [16]. Although 
the square- root variables y 2 and y 3 remain nontrivial in the MRK limit, all of the square roots 
can be rationalized by a clever choice of variables [38]. We define w and w* by 

X = I y = ^ (1 7) 

(i+w)(i+w*y (i + w)(i+w*)' K ' J 

Then the MRK limit of the other variables is 

i 1 1 + w * ~ (l + w)w* 

Ui^-1, 1/1 — >■ 1, V2^y2 = - ri , 2/3^^3 = ^7- tt- (1-8) 

1 + W W{1 + W*) 
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Neglecting terms that vanish like powers of (1 — ui), we expand the remainder function in the 
multi-Regge limit in terms of coefficients multiplying powers of the large logarithm log(l — U\) 
at each loop order, following the conventions of ref. |28j . 




#6(«i,«2,« 3 )|mrk = 2mJ2J2 at l°g n (l-uO [g®(w,w*) + 2mhW{w,w*j\ , (1.9) 



where the coupling constant for planar M = 4 super- Yang-Mills theory is a = g 2 N c / (8tt 2 ). 

The remainder function R$ is a transcendental function with weight 2£ at loop order £. 
Therefore the coefficient functions g$ and hn^ have weight 2£ — n — 1 and 21 — n — 2 respectively. 
As a consequence of eqs. ( 11. 71) and ( 11. 81) . their symbols have only four possible entries, 



Furthermore, w and w* are independent complex variables. Hence the problem of determining 
the coefficient functions factorizes into that of determining functions of w whose symbol entries 
are drawn from {w, 1 + w} — a special class of HPLs — and the complex conjugate functions of 
w*. 

On the other hand, not every combination of HPLs in w and HPLs in w* will appear. When 
the symbol is expressed in terms of the original variables {x, y, y2, 2/3}, the first entry must be 
either x or y, reflecting the branch-cut behavior and first-entry condition for general kinematics. 
Also, the full function must be a single- valued function of x and y, or equivalently a single- valued 
function of w and w*. These conditions imply that the coefficient functions belong to the class 
of SVHPLs defined by Brown [M]. 

The MRK limit (jl.9p is organized hierarchically into the leading-logarithmic approximation 
(LLA) with n = £ — 1, the next-to- leading- logarithmic approximation (NLLA) with n = £ — 2, 
and in general the N fc LL terms with n = £ — k — 1. Just as the problem of DGLAP evolution in 
x space is diagonalized by transforming to the space of Mellin moments N, the MRK limit can 
be diagonalized by performing a Fourier-Mellin transform from (w,w*) to a new space labeled 
by (u, n). In fact, Fadin, Lipatov and Prygarin [38], HO] have given an all-loop-order formula for 
Rq in the multi-Regge limit, in terms of two functions of (z/, n): The eigenvalue u(u,n) of the 
BFKL kernel in the adjoint representation, and the (regularized) MHV impact factor $R eg (i / , n). 
Each function can be expanded in a, and each successive order in a corresponds to increasing k 
by one in the N fc LLA. The leading term in the impact factor is just one, while the leading BFKL 
eigenvalue E v%n was found in ref. [35]. The NLL term in the impact factor was found in ref. |38j . 
and the NLL contribution to the BFKL eigenvalue in ref. [4"0] . 

With this information it is possible to compute the LLA functions gfl^ NLLA functions gf} 2 

and hf_ 2 , and even the real part at NNLLA, hf_ 3 . All one needs to do is perform the inverse 
Fourier-Mellin transform back to the (w,w*) variables. At the three-loop level, this was carried 
out at LLA for # and h? in ref. [38], and at NLLA for g® and in ref. M- Here we 
will use the SVHPL basis to make this step very simple. The inverse transform contains an 
explicit sum over n, and an integral over v which can be evaluated via residues in terms of a sum 




00 e-i 



1=2 n=0 



{w, 1 + W, W*, 1 + W*} . 



(1.10) 
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over a second integer m. For low loop orders we can perform the double sum analytically using 
harmonic sums H7J HBJ HHJ EDI EI] • For high loop orders, it is more efficient to simply truncate 
the double sum. In the (w, w*) plane this truncation corresponds to truncating the power series 
expansion in \ w\ around the origin. We know the answer is a linear combination of a finite number 
of SVHPLs with rational-number coefficients. In order to determine the coefficients, we simply 
compute the power series expansion of the generic linear combination of SVHPLs and match it 
against the truncated double sum over m and n. We can now perform the inverse Fourier-Mellin 
transform, in principle to all orders, and in practice through weight 10, corresponding to 10 loops 
for LLA and 9 loops for NLLA. 

Furthermore, we can bring in additional information at fixed loop order, in order to obtain 
more terms in the expansion of the BFKL eigenvalue and the MHV impact factor. In ref. |40j . 
the NLLA results for g\ and h confirmed a previous prediction [28J based on an analysis of the 

(3) 

multi-Regge limit of the symbol for Rq . In this limit, the two free symbol parameters mentioned 
above dropped out. The symbol could be integrated back up into a function, but a few more 
"beyond-the-symbol" constants entered at this stage. One of the constants was fixed in ref. [10] 
using the NLLA information. As noted in ref. [10], the result from ref. [2H] for can be used 
to determine the NNLLA term in the impact factor. In this paper, we will use our knowledge of 
the space of functions of (w, w*) (the SVHPLs) to build up a dictionary of the functions of (z/, n) 
(special types of harmonic sums) that are the Fourier-Mellin transforms of the SVHPLs. From 

( 31 

this dictionary and we will determine the NNLLA term in the impact factor. 

We can go further if we know the four- loop remainder function . In separate work [52] , we 
have heavily constrained the symbol of R^\ui,u 2 ,u 3 ) for generic kinematics, using exactly the 
same constraints used in ref. [28]: integrability of the symbol, branch-cut behavior, symmetries, 
the final-entry condition, vanishing of collinear limits, and the OPE constraints (which at four 
loops are a constraint on the triple discontinuity). Although there are millions of possible terms 
before applying these constraints, afterwards the symbol contains just 113 free constants (112 if 
we apply the overall normalization for the OPE constraints). Next we construct the multi-Regge 
limit of this symbol, and apply all the information we have about this limit: 

• Vanishing of the super-LLA terms g^ and h^ 1 for n — 4, 5, 6, 7; 

• LLA and NLLA predictions for g^ and hn^ for n = 2, 3; 

• the NNLLA real part h\ , which is also predicted by the NLLA formula; 

• a consistency condition between g^ and . 

Remarkably, these conditions determine all but one of the symbol-level parameters in the MRK 
limit. (The one remaining free parameter seems highly likely to vanish, given the complicated 
way it enters various formulae, but we have not yet proven that to be the case.) 

We then extract the remaining four-loop coefficient functions, g^, and g^\ introducing 
some additional beyond-the-symbol parameters at this stage. We use this information to deter- 
mine the NNLLA BFKL eigenvalue and the N 3 LLA MHV impact factor, up to these parameters. 
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Although our general dictionary of functions of (is, n) contains various multiple harmonic sums, 
we find that the key functions entering the multi-Regge limit can all be expressed just in terms 
of certain rational combinations of v and n, together with the polygamma functions ip, ip', ip", 
etc. (derivatives of the logarithm of the T function) with arguments 1 ± iis + \n\/2. 

As a byproduct, we find that the SVHPLs also describe the multi-Regge limit of the one re- 
maining helicity configuration for six-gluon scattering in M = 4 super- Yang-Mills theory, namely 
the next-to-MHV (NMHV) configuration with three negative and three positive gluon helicities. 
It was shown recently [13] that in LLA the NMHV and MHV remainder functions are related by a 
simple integro-differential operator. This operator has a natural action in terms of the SVHPLs, 
allowing us to easily extend the NMHV LLA results of ref. [43] from three loops to 10 loops. 

This article is organized as follows. In Section [2] we review the structure of the six-point 
MHV remainder function in the multi-Regge limit. Section [3] reviews Brown's construction of 
single- valued harmonic polylogarithms. In Section [4] we exploit the SVHPL basis to determine 
the functions gn and hn at LLA through 10 loops and at NLLA through 9 loops. Section [5J 
determines the NMHV remainder function at LLA through 10 loops. In Section [6] we describe 
our construction of the functions of (is, n) that are the Fourier-Mellin transforms of the SVHPLs. 
Section [7] applies this knowledge, plus information from the four- loop remainder function [52], in 
order to determine the NNLLA MHV impact factor and BFKL eigenvalue, and the N 3 LLA MHV 
impact factor, in terms of a handful of (mostly) beyond-the-symbol constants. In Section [8] we 
report our conclusions and discuss directions for future research. 

We include two appendices. Appendix |A] collects expressions for the SVHPLs (after diago- 
nalizing the action of a Z 2 x Z 2 symmetry), in terms of HPLs through weight 5. It also gives 
expressions before diagonalizing one of the two Z 2 factors. Appendix [B] gives a basis for the func- 
tion space in (is, n) through weight 5, together with the Fourier-Mellin map to the SVHPLs. In 
addition, for the lengthier formulae, we provide separate computer-readable text files as ancillary 
material. In particular, we include files (in Mathematica format) that contain the expressions 
for the SVHPLs in terms of ordinary HPLs up to weight six, decomposed into an eigenbasis of 
the Z 2 x Z 2 symmetry, as well as the analytic results up to weight ten for the imaginary parts 
of the MHV remainder function at LLA and NLLA and for the NMHV remainder function at 
LLA. Furthermore, we include the expressions for the NNLL BFKL eigenvalue and impact factor 
and the N 3 LL impact factor in terms of the building blocks in the variables (is, n) constructed in 
Section [6], as well as a dictionary between these building blocks and the SVHPLs up to weight 
five. 

2 The six-point remainder function in the multi-Regge 
limit 

The principal aim of this paper is to study the six-point MHV amplitude in M = 4 super 
Yang-Mills theory in multi-Regge kinematics. This limit is defined by the hierarchy of scales, 

Sl2 > S345, S456 > S34, S45 , S 56 > S 23 , S 6 i, S 234 . (2.1) 
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In this limit the cross ratios ( II. 2\i behave as 



1 - ux, u 2 , u z ~ , (2.2) 
together with the constraint that the following ratios are held fixed, 

x = = 0(1) and y = = 0{\) . (2.3) 

1 — U\ 1 — Ui 

In the following it will be convenient [38] to parametrize the dependence on x and y by a single 
complex variable w, 

1 ww* . . 

a; = — r and y = — . (2.4) 

(l + w)(l + w*) (l + w)(l + w*) v ; 

Any function of the three cross ratios can then develop large logarithms log(l — «i) in the multi- 
Regge limit, and we can write generically, 

F( Ul , u 2 , u 3 ) = W(l - u i) fi( w > w l + °( l ~ u i) ■ ( 2 - 5 ) 

i 

Let us make at this point an important observation which will be a recurrent theme in the 
rest of the paper: If F(ui,u 2 ,u 3 ) represents a physical quantity like a scattering amplitude, 
then F should only have cuts in physical channels, corresponding to branch cuts starting at 
points where one of the cross ratios vanishes. Rotation around the origin in the complex w 
plane, i.e. (w,w*) — > (e 2m w,e~ 2m w*), does not correspond to crossing any branch cut. As a 
consequence, the functions fi(w, w*) should not change under this operation. More generally, the 
functions fi(w,w*) must be single-valued in the complex w plane. 

Let us start by reviewing the multi-Regge limit of the MHV remainder function R(ui, u 2 , u 3 ) = 
Rq(ui,U2,us) introduced in eq. (II. ip . It can be shown that, while in the Euclidean region the 
remainder function vanishes in the multi-Regge limit, there is a Mandelstam cut such that we 
obtain a non-zero contribution in MRK after performing the analytic continuation [I~6] 

m ^ e~ 2ni \ui\ . (2.6) 

After this analytic continuation, the six-point remainder function can be expanded into the form 
given in eq. (II. 9p . which we repeat here for convenience, 

oo l-\ 

#| M RK = 27riX;Zy log n (l- Ul ) [giP(w,w*)+2mh®(w,w*)] . (2.7) 

1=1 n =0 

The functions gn\w,w*) and hn\w,w*) will in the following be referred to as the coefficient 
functions for the logarithmic expansion in the MRK limit. The imaginary part gn is associated 
with a single discontinuity, and the real part hn with a double discontinuity, although both 
functions also include information from higher discontinuities, albeit with accompanying explicit 
factors of ir 2 . 
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The coefficient functions are single-valued pure transcendental functions in the complex vari- 
able w, of weight 2£ — n — 1 for g$ and weight 2£ — n — 2 for hn- They are left invariant by a 
Z2 x Z2 symmetry acting via complex conjugation and inversion, 

w -H- w* and (w, w*) -H- {1/w, 1/w*) . (2.8) 

The complex conjugation symmetry arises because the MHV remainder function has a parity 
symmetry, or invariance under A — y — A, which inverts y 2 and 2/3 in eq. (11. 8p . The inversion 
symmetry is a consequence of the fact that the six-point remainder function is a totally symmetric 
function of the three cross ratios Ui, u 2 and u 3 . In particular, exchanging y 2 2/3 is the product of 
conjugation and inversion. The inversion symmetry is sometimes referred to as target-projectile 
symmetry [27]. Finally, the vanishing of the six-point remainder function in the collinear limit 

implies the vanishing of gn(w,w*) and h„ (w,w*) in the limit where (w,w*) — > 0. Clearly the 

(£) (£) 

functions g n and h n are already highly constrained on general grounds. 

In ref. [381 SO] an all-loop integral formula for the six-point amplitude in MRK was presented^, 

00 » f+00 r / 1 \ w[y,n) 

e R ^% KK = C o^uJ ab + ^ a - T, (-IT (-Y I -^^H^W^n) ' 

V V «2 "3 / 

(2.9) 



2 „^L Kw J J-00 v 2 + \ \ y/u^i 



The first term is the Regge pole contribution, with 

Uab = \ ik{o) log — = -7k(o) log I w 1 2 , (2.10) 
8 u 2 8 

and 7_ftr(a) is the cusp anomalous dimension, known to all orders in perturbation theory [53] . 

00 

lK {a) = J2^K a£ = 4a-4C 2 a 2 + 22C4a 3 -(^C6 + 4C 3 2 )a 4 + ---. (2.11) 

1=1 

The second term in eq. (12.91) arises from a Regge cut and is fully determined to all orders by 
the BFKL eigenvalue u(u,n) and the (regularized) impact factor $R eg (r/, n). The function 5 
appearing in the exponent on the left-hand side is the contribution from a Mandelstam cut 
present in the BDS ansatz, and is given to all loop orders by 

1 1 \w\ 2 

S = - 1k{o) log (xy) = - <y K (a) log 4 . (2.12) 



In addition, we have 

1 1 \l + w\ 



y/U 2 U 3 1 - Ui \w\ 



(2.13) 



1 There is a difference in conventions regarding the definition of the remainder function. What we call R is 
called log(i?) in refs. [351 HD]- Apart from the zeroth order term, the first place this makes a difference is at four 
loops, in the real part. 
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The BFKL eigenvalue and the impact factor can be expanded perturbatively, 
u(v, n) = -a (E„ tn + a + a 2 + 0(a 3 )) , 



$ Rcg (z/, n) = 1 + a $£ g (z/, n) + a' Sg^i/, n) + a 3 $^ g (z/, n) + (9(a 4 ) . 
The BFKL eigenvalue is known to the first two orders in perturbation theory (301 ES], 



(2) 



,(3) 



(2.14) 



E,. 



E^i 



\n\ 



+ iff[l + iu + 



\n\ 
~2 



+ V 1 



l"l 

T 



2V(1) 



1 

2ii/ 



21/ 



|n| 

T 



+ ^" l - ii/ + 



|n| 

T 



(2.15) 
(2.16) 



^ ( 1 + iv + y ] - $ I 



2Z/ 



1 71 1 

~2 



, 1 \ n \ 

3C3-T 



" - T 



4 (^ + if) 3 
7s is the Euler-Mascheroni 

(2.17) 



where ip{z) = j-logr(z) is the digamma function, and ip(l) = 
constant. The NLL contribution to the impact factor is given by 



2 v ' n 



n 



y + f) 2 



(2 



The BFKL eigenvalues and impact factor in eqs. (12.151) . (12.161) and (I2.17P are enough to compute 
the six-point remainder function in the Regge limit in the leading and next-to-leading logarithmic 
approximations (LLA and NLLA). Indeed, we can interpret the integral in eq. (12. 9 j) as a contour 
integral in the complex v plane and close the contour at infinity. By summing up the residues we 
then obtain the analytic expression of the remainder function in the LLA and NLLA in MRK. 
This procedure will be discussed in greater detail in Section |H Some comments are in order 
about the integral in eq. (12. 9ft : 

1. The contribution coming from n = seems ill-defined, as the integral in eq. (12. 9 j) diverges. 
After closing the contour at infinity, our prescription is to take only half of the residue at 
v = n = into account. 

2. We need to specify the Riemann sheet of the exponential factor in the right-hand side of 
eq. (I2.9p . We find that the replacement 



y/U 2 U 3 



u)(y,n) 



-iTruj(is,n) 



V M 2 U 3 



(2.18) 



gives the correct result. 

The in factor in the right-hand side of eq. (I2.18P generates the real parts hn in eq. (12. 7p . It 
is easy to see that the g£) and hffl functions are not independent, but they are related. For 
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example, at LLA and NLLA we have, 
hfl^w,™*) =0, 

i *\ (t) f *\ , 1 (i) V-i)/ *\ i |1 + H 4 

h\l 2 [w,w) = —^-g\l x {w,w ) + —Ik 9e-2 K w ) lo S | w | 2 ^ ig ^ 

where 7^ = 4 from eq. (12.111) . (Note that the sum over k in the formula for hfl 2 would not 
have been present if we had used the convention for R in refs. [3S1 SO]-) Similar relations can be 
derived beyond NLLA, i.e. for n < i — 2. 

So far we have only considered 2 — > 4 scattering. In ref. [39] it was shown that if the remainder 
function is analytically continued to the region corresponding to 3 — > 3 scattering, then it takes 
a particularly simple form. The analytic continuation from 2 — > 4 to 3 — > 3 scattering can be 
obtained easily by performing the replacement 

log(l - Mi) log(«i — 1) — *7T (2.20) 

in eq. (I2.9p . After analytic continuation the real part of the remainder function only gets contri- 
butions from the Regge pole and is given by [39] 

Re ( e ^- i7n5 ) = cos nco ab . (2.21) 

It is manifest from eq. (12.91) that eq. (12.211) is automatically satisfied if the relations among the 
coefficient functions derivable by tracking the in from eq. (I2.18P (e.g. eq. (I2.19P ) are satisfied in 
2 — ?■ 4 kinematics. 

So far we have only reviewed some general properties of the six-point remainder function in 
MRK, but we have not yet given explicit analytic expressions for the coefficient functions. The 
two- loop contributions to eq. (12.91) in LLA and NLLA were computed in refs. [371 ESj , while the 
three-loop contributions up to the NNLLA were found in refs. [281 EZ]- 111 all cases the results 
have been expressed as combinations of classical polylogarithms in the complex variable w and 
its complex conjugate w*, with potential branching points at w = and w — — 1. As discussed 
at the beginning of this section, all the branch cuts in the complex w plane must cancel, i.e., the 
function must be single-valued in w. The class of functions satisfying these constraints has been 
studied in full generality in the mathematical literature, as will be reviewed in the next section. 



3 Harmonic polylogarithms and their single-valued ana- 
logues 

3.1 Review of harmonic polylogarithms 

In this section we give a short review of the classical and harmonic polylogarithms, one of the 
main themes in the rest of this paper. The simplest possible polylogarithmic functions are the 



11 



so-called classical polylogarithms, denned inside the unit circle by a convergent power series, 



fc=i 



They can be continued to the cut plane C\[l, oo) by an iterated integral representation, 



U m (z)= I' dz' Um -] {Z,) . (3.2) 
Jo z 

For m — 1, the polylogarithm reduces to the ordinary logarithm, Lii(z) = — log(l — z), a fact 
that dictates the location of the branch cut for all m (along the real axis for z > 1). It also 
determines the discontinuity across the cut, 

log™" 1 z 

ALi m (z) = 2m S . (3.3) 
[m — ly. 

It is possible to define more general classes of polylogarithmic functions by allowing for dif- 
ferent kernels inside the iterated integral in eq. (13. 1 1) . The harmonic polylogarithms (HPLs) jl5] 
are a special class of generalized polylogarithms whose properties and construction we review in 
the remainder of this section. To begin, let w be a word formed from the letters Xq and x±, and 
let e be the empty word. Then, for each w, define a function H w (z) which obeys the differential 
equations, 

—H X0W {z) = ^^- and = , ( ^ 

subject to the following conditions, 

H e {z) = l, H x n{z) = ^log n z, and lim H w ^(z) = . (3.5) 

There is a unique family of solutions to these equations, and it defines the HPLs. Note that we use 
the term "HPL" in a restricted sensed - we only consider poles in the differential equations f l3.4[) 
at z = and z = 1. (In our MRK application, we will let z = —w, so that the poles are at w = 
and w = —1.) 

The weight of an HPL is the length of the word w, and its depth is the number of xi'dE 
HPLs of depth one are simply the classical polylogarithms, H n (z) = Li n (z). Like the classical 
polylogarithms, the HPLs can be written as iterated integrals, 

H xow (z) = / dz — and H X1W = / dz f- . (3.7) 

Jo z Jo 1 — z 



2 In the mathematical literature, these functions are sometimes referred to as multiple polylogarithms in one 
variable. 

3 For ease of notation, we will often impose the replacement {xq — > 0, x\ — > 1} in subscripts. In some cases, we 
will use the collapsed notation where a subscript m denotes m — 1 zeroes followed by a single 1. For example, if 

H w (z) = H XQXoXlXoXl (z) = H 0t 0,1,0,1 (z) = H 3t2 (z) . (3.6) 
In the collapsed notation, the weight is the sum of the indices, and the depth is the number of nonzero indices. 
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Weight 


Lyndon words 


Dimension 


1 


0, 1 


2 


2 


01 


1 


3 


001, Oil 


2 


4 


0001, 0011, 0111 


3 


5 


00001, 00011, 00101, 00111, 01011, 01111 


6 



Table 1: All Lyndon words Lyndon(x , X\) through weight five 



The structure of the underlying iterated integrals endows the HPLs with an important property: 
they form a shuffle algebra. The shuffle relations can be written, 

H wi (z) H W2 (z) = H ^ z ^ ( 3 - 8 ) 

where u>iirru;2 is the set of mergers of the sequences W\ and W2 that preserve their relative 
ordering. Equation ( 13. 8 p may be used to express all HPLs of a given weight in terms of a 
relatively small set of basis functions and products of lower-weight HPLs. One convenient such 
basis [M] of irreducible functions is the Lyndon basis, defined by {H w (z) : w £ Lyndon(x , 
The Lyndon words Lyndon(xo, X\) are those words w such that for every decomposition into two 
words w = uv, the left word is lexicographically smaller than the right, u < v. Table [1] gives the 
first few examples of Lyndon words. 

All HPLs are real whenever the argument z is less than 1, and so, in particular, the HPLs 
are analytic in a neighborhood of z = 0. The Taylor expansion around z = is particularly 
simple and involves only a special class of harmonic numbers [HJ HH] (hence the name harmonic 
polylogarithm), 

00 ; 

H mi ^ mk (z) = ^—Z rn2t ... >mk (l-1), m;>0, (3.9) 
where Z mu _ >mk (n) denote the so-called Euler-Zagier sums (HUE?], defined recursively by 

n -y n 1 

Z mi (n) = ^ — ■ and Z mi ,..., mfe (n) = ^ — Z m2 ^ mk (l - 1) . (3.10) 
i=i 1=1 

Note that the indexing of the weight vectors mi, . . . , in eqs. (13. 9p and (I3.10p is in the collapsed 
notation. 

Another important property of HPLs is that they are closed under certain transformations of 
the arguments [IS]- In particular, using the integral representation (13. 7p . it is easy to show that 
the set of all HPLs is closed under the following transformations, 

z 1 — y 1 — Zj z I—)- 1/z, z \-> 1/(1 - z), z^l — l/z, z^zjiz — V). (3-11) 

If we add to these mappings the identity map z 1— > z, we can identify the transformations in 
eq. ( 13. lip as forming a representation of the symmetric group S3. In other words, the vector 
space spanned by all HPLs is endowed with a natural action of the symmetric group S3. 
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Finally, it is evident from the iterated integral representation (13. 7p that HPLs can have branch 
cuts starting at z = and/or z = 1, i.e., HPLs define in general multi- valued functions on the 
complex plane. In the next section we will define analogues of HPLs without any branch cuts, 
thus obtaining a single-valued version of the HPLs. 



3.2 Single- valued harmonic poly logarithms 

Before reviewing the definition of single-valued harmonic polylogarithms in general, let us first 
review the special case of single- valued classical polylogarithms. The knowledge of the disconti- 
nuities of the classical polylogarithms, eq. (13. 3p . can be leveraged to construct a sequence of real 
analytic functions on the punctured plane C\{0, 1}. The idea is to consider linear combinations 
of (products of) classical polylogarithms and ordinary logarithms such that all the branch cuts 
cancel. Although the space of single- valued functions is unique, the choice of basis is not unique, 
and there have been several versions proposed in the literature. As an illustration, consider the 
functions of Zagier [55J, 



k=l 



where 9i m denotes the imaginary part for m even and the real part for m odd. The discontinuity 
of the function inside the curly brackets is given by 

( -;° gi irif~ ,\i = v^w**')- 1 - < 3 - i3 > 

^ (m — k)\ (A; — 1)! (m — 1)! 



k=l 



Since eq. ( 13 . 1 3[) is real for even m and pure imaginary for odd m, D m (z) is indeed single- valued. 
For the special case m = 2, we reproduce the famous Bloch-Wigner dilogarithm [56] . 



D 2 (z) = Im{Li 2 (^)} + arg(l - z) log \z\ . (3.14) 

Just as there have been numerous proposals in the literature for single-valued versions of the 
classical polylogarithms, there are many potential choices of bases for single-valued HPLs. On 
the other hand, if we choose to demand some reasonable properties, it turns out that a unique set 
of functions emerges. Following ref. [44], we require the single- valued HPLs to be built entirely 
from holomorphic and anti-holomorphic HPLs. Specifically, they should be a linear combination 
of terms of the form H wl (z)H W2 (z), where W\ and w 2 are words in xo and x\ or the empty 
word e. The single-valued classical polylogarithms obey an analogous property, and it can be 
understood as the condition that the single-valued functions are the proper extensions of the 
original functions. The remaining requirements are simply the analogs of the conditions used to 
construct the ordinary HPLs. 

Define a function C w (z), which is a linear combination of functions H wl (z)H W2 (z) and which 
obeys the differential equations 

d £w(z) ^ C w iz) 

fcZxowW = and ~Q~z XlW ^ = T^z~' 
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subject to the conditions, 

£ e 0) = l, C x n(z) = ^\og n \z\ 2 and \im£ w ^ x n(z) = . (3.16) 

In ref. jH] Brown showed that there is a unique family of solutions to these equations that is 
single- valued in the complex z plane, and it defines the single- valued HPLs (SVHPLs). The 
functions C w (z) are linearly independent and span the space. That is to say, every single- valued 
linear combination of functions of the form H Wl (z)H W2 (z) can be written in terms of the C w (z). 
In ref. [H] an algorithm was presented that allows for the explicit construction of all SVHPLs as 
linear combinations of (products of) ordinary HPLs. We present a short review of this algorithm 
in Section I3TB1 

The SVHPLs of ref. [2] share all the nice features of their multi- valued analogues. First, like 
the ordinary HPLs, they obey shuffle relations, 

c wi (z) c W2 (z) = ^2 £w(*0> ( 3 - 17 ) 

where again wiiiiw^ represents the shuffles of W\ and wi. As a consequence, we may again choose 
to solve eq. (I3.17P in terms of a Lyndon basis. It follows that if we want the full list of all SVHPLs 
of a given weight, it is enough to know the corresponding Lyndon basis up to that weight. 

Furthermore, the space of SVHPLs is also closed under the S3 action defined by eq. (13. lip . 
Indeed, if we extend the action to the complex conjugate variable z, then the closure of the 
space of all ordinary HPLs implies the closure of the space spanned by all products of the form 
H wl (z)H W2 (z), and, in particular, the closure of the subspace of SVHPLs. For the SVHPLs, it 
is possible to enlarge the symmetry group to Z 2 x S 3 , where the Z 2 subgroup acts by complex 
conjugation, z -H- z. 

It turns out that the functions C w (z) can generically be decomposed as 

C w (z) = (H w (z) - (-l) lwl H w {z)) + [products of lower weight] , (3.18) 

where \w\ denotes the weight. As such, it is convenient to consider the even and odd projections, 
i.e., the decomposition into eigenfunctions of the Z 2 action, 

L w (z) = l - (£ w (z) - (-1)^1 C w (z)) and L w (z) = ~ (C w (z) + (-lpC w (z)) . (3.19) 

The basis defined by C w (z) was already complete, and yet here we have doubled the number of 
potential basis functions. Therefore L w (z) and L w (z) must be related to one another. Writing 
L w (z) = 9l\ w \(£ w (z)) , we see that it has the same parity as Zagier's single-valued versions of 
the classical polylogarithms given in eq. f)3.12p . Therefore we might expect the L w (z) to form a 
complete basis on their own. Indeed this turns out to be the case, and the L w (z) can be expressed 
as products of the functions L w (z), 

L w (z) = [products of lower weight L w r(z)] . (3.20) 
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Hence we will not consider the functions L w (z) any further and will concentrate solely on the 
functions L w (z). 

The functions L w (z) do not automatically form simple representations of the S3 symmetry. 
For the current application, we will mostly be concerned with the Z 2 C S3 subgroup generated 
by inversions z 1/z. The functions L w (z) can easily be decomposed into eigenfunctions of this 
Z 2 , and, furthermore, these eigenfunctions form a basis for the space of all SVHPLs. The latter 
follows from the observation that, 

L w (z) - (-1)M+*»£, W (_) = [products of lower weight], (3.21) 

where \w\ is the weight and d w is the depth of the word w. We will denote these eigenfunctions 
of Z 2 x Z 2 by, 

1 



L 



w 



'z) ± L w Q-_ 



(3.22) 



and present most of our results in terms of this convenient basis. For low weights, Appendix |A] 
gives explicit representations of these basis functions in terms of HPLs. The expressions through 
weight six can be found in the ancillary files. 

We have seen in the previous section that in the multi-Regge limit the six-point amplitude is 
described to all loop orders by single-valued functions of a single complex variable w satisfying 
certain reality and inversion properties. It turns out that the SVHPLs we just defined are 
particularly well-suited to describe these multi-Regge limits. This description will be the topic 
of the rest of this paper. 



3.3 Explicit construction 

The explicit construction of the functions C w (z) is somewhat involved so we take a brief detour 
to describe the details. Let X* be the set of words in the alphabet {xq, along with the empty 
word e. Define the Drinfel'd associator Z(xq, x\) as the generating series, 

Z(x ,x 1 ) = ^ C(w)w, (3.23) 
wex* 

where ((w) = H W {1) for w 7^ x± and C( x i) — 0. Using the collapsed notation for w, these ((w) 
are the familiar multiple zeta values. 

Next, define an alphabet {yo,yi} (and a set of words Y*) and a map ~ : Y* — > Y* as the 
operation that reverses words. The alphabet {yo,yi} is related to the alphabet {xo,xi} by the 
following relations: 

yo = xq 

1 _i (3-24) 

Z(yo,yi)yiZ{y ,yi) = Z(x ,xi) xiZ(x 0l xi). 

The inversion operator is to be understood as a formal series expansion in the weight \w\. Solving 
eq. ( I3.24p iteratively in the weight yields a series expansion for yi. The first few terms are, 

!/1 = I1 -C 3 (2 WlIl -4 IoWl + 2 Wia;i+ 4 Wia;o 

— 6xiX XiXi — 2xiXiX x + 6xiXiX Xi — 2xio:ia;ia;o) + . . . 
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Letting : Y* — >■ X* be the map that renames y to x, i.e. (j)(yo) = %o and = X\, define the 
generating functions 

Lx(z) = ^2 H w(z)w, L Y (z) = ^ H^ w )(z)w. (3.26) 

weX* w<=Y* 

In the following, we use a condensed notation for the HPL arguments, in order to improve the 
readability of explicit formulas: 

H w = H w (z) and H w = H w {z) . (3.27) 

Then we can write 

L x (z) = 1 + H xq + Hxxx 



+ H 0fi x x + Hq^xqXi + H lfi xix + H 1;1 xixi 

+ -^0,0,0^0^0^0 + -^0,0,1^0^0^1 + -^0,1,0^0^1^0 + -f^O.ljl^O^l^l 

+ HiflflXix^XQ + Hia^xiXqXi + i/i 5 i i0 a;iXiXo + #1,1,1X1X1X1 + 



(3.28) 



and 



(3.29) 



L Y (z) = 1 + H y + H l y l 

+ H 0fi y y + H ,iyiy + -ff~i,o2/o2/i + #i,i2/i2/i 
+ H 0fifi y y y + # ,o,i2/i2/o?/o + #o,i,oZ/oZ/iyo + #o,i,i2/i2/i2/o 
+ H lfifi y y y 1 + H 1A1 y 1 y y 1 + H 1 , lfi y y 1 y 1 + H\^\y\y\yi + ... 

= 1 + #0^0 + -fflXi 

+ i7 ,oXox + -f/o.i^o + Hi.o^o^i + -f^l,lXiXi 

+ ^0,0,0^0^0^0 + ^0,0,1^1^0^0 + ^0,1,0^0^1^0 + ^0,1,1^1^1^0 

+ #1,0,0^0^0^1 + -^1,0,1^1^0^1 + -H" 1,1,0^0^1^1 + ^1,1,1^1^1^1 + • • • • 

In the last step of eq. ( I3.29|) we used yo = xq and yi — x±. Note that the latter only holds 
through weight three, as is clear from eq. (I3.25p . Finally, we are able to construct the SVHPLs 
as a generating series, 

C(z) = L x (z)L Y (z) = C w (z)w. (3.30) 

w&X* 

Indeed, taking the product of eq. (I3.28P with eq. f!3.29j) and keeping terms through weight three, 
we obtain, 

s ^C w (z)w = 1+ £ (z) x + £i(z) xi 

wex* 

+ £0,0(2) x x + £0,1(2) x xi + £1,0(2) xix + £1,1(2) XiX! (3.31) 

+ £o,0,o(^) X0X0X0 + £0,0,1(2) X0X0X1 + £0,1,0(2) XoX\Xo + £0,1,1(2;) X0X1X1 

+ £1,0,0(2) XiX X + £l,0,l(2) XiX Xi + £1,1,0(2) X1X1X0 + £1,1,1(2) £1X1X1 + . . . , 
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where the SVHPL's of weight one are, 

C (z) = Ho + Ho, Ci(z) = Hx + Hx, 
the SVHPL's of weight two are, 

£0,0(2) = #0,0 + #0,0 + H H , 

£0,1(2) = #0,1 + H lt o + H Hx , 

£i,o(*) = #i,o + #o,i + #i#o, 

Cx,x{z) = Hx A + Hx,x + HxHx, 
and the SVHPL's of weight three are, 



(3.32) 



(3.33) 



£0,0,0(2; 
£0,0,1(2; 

£0,1,0(2; 
£0,1,1(2; 
£1,0,0(2; 
£1,0,1(2; 
£1,1,0(2; 
£1,1,1(2; 



#0,0,0 + -"0,0,0 + #o,o#o + HoHofi , 
#0,0,1 + #1,0,0 + #o,o#i + H Hx,o , 
Ho,x,o + H ,x,o + H 0: iH + H H 0;1 , 
-f^o,i,i + ^1,1,0 + H 0>1 Hx + H H lt x , 
#1,0,0 + #0,0,1 + #i,o#o + #i#o,o , 
#1,0,1 + #1,0,1 + #i,o#i + #i#i,o 5 
#1,1,0 + #0,1,1 + #i,i#o + #i#o,i , 
#1,1,1 + #1,1,1 + #i,i#i + #i#i,i • 



(3.34) 



The 1/ alphabet differs from the x alphabet starting at weight four. Referring to eq. (I3.25p . we 
expect the difference to generate factors of (3. To illustrate this effect, we list here the subset of 
weight-four SVHPLs with explicit ( terms: 



£0,0,1,1(2; 
£0,1,0,1(2; 
£0,1,1,1(2; 
£1,0,1,0(2; 
£1,0,1,1(2; 
£1,1,0,0(2 
£1,1,0,1(2 
£1,1,1,0(2; 



#0,0,1,1 + #1,1,0,0 + #0,0,1#1 + #0#1,1,0 + #o,o#i,i — 2C3 #1 

#0,1,0,1 + #1,0,1,0 + #0,1, 0#1 + #0#1,0,1 + #0,1#1,0 + 4C3 #1 

#0,1,1,1 + #1,1,1,0 + #0,1, 1#1 + #0#1,1,1 + #o,i#i,i — 2C3 #1 

#1,0,1,0 + #0,1,0,1 + #1,0,1#0 + #1 #0,1,0 + #i,o#o,i — 4C3 #1 

#1,0,1,1 + #1,1,0,1 + #1,0,1#1 + #1 #1,1,0 + #1,0#1,1 + 6C3 #1 

#1 ,1,0,0 + #0,0,1,1 + #1,1, 0#0 + #1 #0,0,1 + #1,1#0,0 + 2(^3 #1 

#1,1,0,1 + #1,0,1,1 + #1,1, 0#1 + #1#1,0,1 + #i,i#i,o — 6C3 #1 

#1,1,1,0 + #0,1,1,1 + #1,1, 1#0 + #1 #0,1,1 + #1,1#0,1 + 2C3 #1 



(3.35) 



Finally, we remark that the generating series £(2;) provides a convenient way to represent the 
differential equations (I3.15p . Together with the y alphabet, it also allows us to write down the 
differential equations in z, 



d_ 

dz 



C(z) 



z 



Xx 

l-z 



C{z) 



and 







dz 



-C{z) = C(z) 



yo 

z 



IJl 



(3.36) 



These equations will be particularly useful in Section when we study the multi-Regge limit of 
the ratio function of the six-point NMHV amplitude. 



4 The six-point remainder function in LLA and NLLA 

In Section [21 we showed that in MRK the remainder function is fully determined by the coefficient 
functions gn(w,w*) and hn\w,w*) in the logarithmic expansion of its real and imaginary part 
in eq. (12 .7p . We further argued that these functions are single- valued in the complex w plane, 
and suggested that they can be computed explicitly by interpreting the z/-integral in eq. (12. 9p as 
a contour integral and summing the residues. In this section, we describe how knowledge about 
the space of SVHPLs can be used to facilitate this calculation. In particular, we present results 
for LLA through ten loops and for NLLA through nine loops. 

The main integral we consider is eq. (12. 9p . which we reproduce here for clarity, rewriting the 
last factor to take into account eqs. (12. 13[) and (I2.18p . 



OO n /* -|-(X) 7 

e* + -W = + 4 £ (-1)" (^Y / -^M 2w W^ 



x exp 



—uj(u, n) log(l - ui) + in + - lo 



1 , \w 



2 



2 °|l + w| 4 



• (4-1) 



The integrand depends on the BFKL eigenvalue and impact factor, which are known through 
order a 2 and are given in eqs. (12. 15ft . (12. 16ft and (12.171) . These functions can be written as 
rational functions of v and n, and polygamma functions (ip and its derivatives) with arguments 
\n\/2. Recalling that the polygamma functions have poles at the non-positive integers, it 
is easy to see that all poles are found in the complex v plane at values v = — i{m + 4^), m G N, 
n G Z. When the integral is performed by summing residues, the result will be of the form, 



E 



a m , n w m+ "-w' m . (4.2) 



Because residues of the polygamma functions are rational numbers, and because polygamma 
functions evaluate to Euler-Zagier sums for positive integers, the coefficients a m ^ n are combina- 
tions of 

1. rational functions in m and n, 

2. Euler-Zagier sums of the form Z$(m), Z?(n) and Z^im + n), 

3. log \ w\, arising from taking residues at multiple poles. 

Identifying (z, z) = (—w, —w*), and comparing the double sum (14. 2 p to the formal series expan- 
sion of the HPLs around z = 0, eq. (13. 9p . we conclude that the double sums will evaluate to 
linear combinations of terms of the form H Wl (w)H W2 (w*). Moreover, as discussed above, this 



19 



combination should be single- valued. Therefore, based on the discussion in Section [21 we expect 



and hn\w,w*) to belong to the space spanned by the SVHPLs. 



Furthermore, we know that g4 (w,w*) and hn'(w,w*) are invariant under the action of the 
Z 2 x Z 2 transformations of eq. ( 12.81) . In terms of SVHPLs, this symmetry is just an (abelian) 
subgroup of the larger Z 2 x S3 symmetry, where the Z 2 is complex conjugation and the S3 action 
is given in eq. (13. lip . As such, we do not expect an arbitrary linear combination of SVHPLs, 
but only those that are eigenf unctions with eigenvalue (+, +) of the Z 2 x Z 2 symmetry. 

Putting everything together, and taking into account that scattering amplitudes in M = 4 
SYM are expected to have uniform transcendentality, we are led to conjecture that, to all loop 
orders, gn\w,w*) and hn (w,w*) should be expressible as a linear combination of SVHPLs in 
(z, z) = (—w, —w*) of uniform transcendental weight, with eigenvalue (+, +) under the Z 2 x Z 2 
symmetry. Inspecting eq. ( 12. 7|) , the weight should be 2£ — n — 1 for g n and 2£ — n — 2 for h n . 
Our conjecture allows us to predict a priori the set of functions that can appear at a given loop 
order, and in practice this set turns out to be rather small. Knowledge of this set of functions 
can be used to facilitate the evaluation of eq. (14. ip . We outline two strategies to achieve this: 

1. Evaluate the double sum (14.21) with the summation algorithms of ref. [57]. The result is 
a complicated expression involving multiple polylogarithms which can be matched to a 
combination of SVHPLs and zeta values by means of the symbol [20| I2TI [22| [23| 124"] and 
coproduct [SHUSHED]. 



2. The double sum (14 .2p should be equal to the formal series expansion of some linear combi- 
nation of SVHPLs and zeta values. The unknown coefficients of this combination can be 
fixed by matching the two expressions term by term. 



To see how this works, we calculate the two-loop remainder function in MRK. Expanding eq. ( 14. ip 
to two loops, we find, 



■-L+ + -Z[l] 
2 1 4 



+ 



log(l - Ul ) \x[E v , n ] + ^-C 2 Lt + -X^M] + \l\X\E v , 
+2m (1 [L ] 2 + ~ [Lf] 2 - I LtX[l] + \x[E v , n ] 



(4.3) 



where we have introduced the notation, 



7T *■ — ' 



w \ 2 / dv 



IV* 



2 1 VL 

00 ' A 



(4.4) 



Explicit expressions for the functions for low weights are provided in Appendix |A] Equa- 
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tion (14. 3 p is consistent only if the term of order a vanishes. Indeed this is the case, 

s r+oo dv 



4 



2iu 



oo V + — 



= , ogM 2 + 2 v^: + 2f: tmT (4.5) 

n=l n=l 

= log \w\ 2 — 2 log 1 1 + w\ 2 
= 2Lf . 

As previously mentioned, we only take half of the residue at v — n = 0. 

Moving on to the terms of order a 2 , we refer to eq. (12.71) and extract from eq. (14.31) the 
expressions for the coefficient functions, 

gi 2 \w,w*) = ±Z[E v , n ] 

gi 2 \w,w*) = \<M + \l[$ZM] + \ L tl[Ev,n] (4.6) 
h?\w,w*) = 1 [L -] 2 + 1 [L+] 2 - Il+J[1] + i J[J^ lB ] . 

Note that /i^ = 0, in accordance with the general expectation that hf^ = 0. Proceeding 
onwards, we have to calculate X[E Vtn ], 



*J=L Kw J J '-00 ^ 2 + ! r I 



2(z/ 2 + 



-,.2 - 



71 \ / 71 



- r h£ _ 2 (-u.)- + (-u.-)-' + p |w|J + (m)] (-«,)- + (-■»' 

m=l v 



77 i 



(4.7) 



m m + n 

n=l m=l v ' 

The single sum in the first line immediately evaluates to polylogarithms, 

H 2m J-w) m + (-w*) m 2 (-w) m + (-w*y 



E(^- 2 '-"'f''' + Hog M 2 + 2ft W] 
^— ' I m 2 ttt 

m=l v 



m 

V + log w 2 + 2Z 1 (m - 1 P ^ ^ '- 

m=i ' 

log |uf [ifi(-w) + Hx{-w*)} + 2# ,i(M 2 ) + 2H ltl {-w) + 2H ljl {-w*) . 
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Next we transform the double sum into a nested sum by shifting the summation variables by 
n = N — m, 



n=l m=l v y AT=1 m=l < 

= Li^^-w, —it;*) + Lii jl (— —w) 

_ U I ZJ I „.,*\ U / 1 „ .. 1 2\ 



(4.9) 



where the last step follows from a stuffle identity among multiple polylogarithms [HI]. Putting 
everything together, we obtain 



X[E u>n ] = log \w\ 2 [Hxi-w) + Hti-w*)] + 2H ltl {-w) + 2H^{-w*) + 2H 1 (-w) Hx(-w*) 

(4.10) 



Referring to eqs. H4.5[) and H4.6[) , we can now write down the results, 

yi v > ; 4 16 1 0J (4.11) 

h { o\w,w*) =0. 

For higher weights the nested double sums can be more complicated, but they are always 
of a form that can be performed using the algorithms of ref. [57]. These algorithms will in 
general produce complicated multiple polylogarithms that, unlike in eq. ( I4.9p . cannot in general 
be reduced to HPLs by the simple application of stuffle identities. In this case we can use 
symbols [221 1231 121] and the coproduct on multiple polylogarithms [SHI ESI ED] to perform this 
reduction. 

The above strategy becomes computationally taxing for high weights. For this reason, we also 
employ an alternative strategy, based on matching series expansions, which is computationally 

(2) 

simpler. We demonstrate this method in the computation of g$ , for which the only missing 
ingredient in eq. (14. 6 p is Z[$j^ n)], where 3>Reg( z/ > n ) * s defined in eq. (I2.17p . To proceed, we 
write the //-integral as a sum of residues, and truncate the resulting double sum to some finite 
order, 




(4.12) 
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< 2 log \w\ 2 - (log |u>| 2 ) 



w\ 



1 + - log \w\ 
4 



+ (w + iw*) 

i / 2 i *2\ 

+ . . . . 



2C2 



4 — 2 log | | 2 H — log \w\ 



\w\* + ... 



1 



1 + ^log \w\ 2 



M + ... 



1 X l 2, 

- + -log |w| 



-1 log Iwl 2 

3 &l 1 



Here we show on separate lines the contributions to the sum from n = 0, n = ±1, and n = ±2. 
Next, we construct an ansatz of SVHPLs whose series expansion we attempt to match to the 
above expression. We expect the result to be a weight-three SVHPL with parity (+, +) under 
conjugation and inversion. Including zeta values, there are five functions satisfying these criteria, 
and we can write the ansatz as, 



Z[$S>, n)] = ci Lt + c 2 [L ] 2 Lt + c 3 [L+f + c 4 C2 ^ + c 5 Cs • 



-l2r + 



(4.13) 



Using the series expansions of the constituent HPLs ( 13.9f) . it is straightforward to produce the 
series expansion of this ansatz, 



c l c 2 c 3 

12 2" 8" 



log 3 |w;| 2 + -C4C2 log M 2 + c 5 C 3 + 3 c 3 (log |w| 2 ) L 



10 



+ (w + w*) 

+ . . . . 



Cl 



3c 3 



C2C4 + -Ci + -ci log |w| + — - - c 2 — log \w\ + ... 



4 



(4.14) 



We have only listed the terms necessary to fix the undetermined constants. In practice we 
generate many more terms than necessary to cross-check the result. Consistency of eqs. (I4.12p 
and (I4.14p requires, 



Cl 



c 2 



c 3 



1 

3' 



-2, c 5 = 



which gives, 



^>£ g M] = -4^ 3 + + \ [L ] 2 Lt - I [L+f 



Finally, putting everything together in eq. ( I4.6j) . 



(2)/ * 
g y '(w,w 



-Lt + l[Lff + i[L ] 2 L+ 



(4.15) 



(4.16) 



(4.17) 



This completes the two- loop calculation, and we find agreement with (33, [38]. Moving on to 
three loops, we can proceed in exactly the same way, and we reproduce the LLA [38] and NLLA 
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results J2SJ HO] for the imaginary parts of the coefficient functions, 

g{ 3 \w, w*) = - §tf Li + ^[Ltr + f [Lt? - ^[Lo]' (4-18) 

IT 2 IT 2 1 

12 L 1J 4g L OJ T 4 Sd i 

(The result for gj 3 '' agrees with that in ref. once the constants are fixed to c = and 
7' = —9/2 [30].) The real parts are given by, 

f$\w,w*) =0, 

™ 111 ( 4 - 19 ) 

in agreement with ref. [3B]- Using the fact that 

1 \w\ 2 

L+ = - log , 11 j . , (4.20) 
1 2 & |l + u;| 4 ' V ; 

it is easy to check that hf\w, w*) satisfies eq. (I2.19P for £ = 3. 

It is straightforward to extend these methods to higher loops. We have produced results for 
all functions with weight less than or equal to 10, which is equivalent to 10 loops in the LLA, 
and 9 loops in the NLLA. Using the C++ symbolic computation framework GiNaC [62], which 
allows for the efficient numerical evaluation of HPLs to high precision [63] , we can evaluate these 
functions numerically. Figures [1] and [2] show the functions plotted on the line segment for which 
w = w* and < w < 1. Here we also show the analytical results through six loops. We provide 
a separate computer-readable text file, compatible with the Mathematica package HPL [611 EH], 
which contains all the expressions through weight 10. 

Up to six loops, we find, 

= ^[L2} 2 + ^[Lo] 2 [Li} 2 + ^[L^ + ^[Ltr-^L^L^ (4.21) 



= ^[Lo] 2 [Lt] 3 + ^Lt[L^-^Li[L ] 2 + 1 -[L^] 2 C 3 (4.22) 



48 1 2 J 48 




-48^ L3+ " 


8 lCs ' 


l^o mi 3 


+ i4 Ll 


-^] 2 Cs + 




3 r+ 1 




^8 3,111 4 


L+ 5 



24 L U 16 2>1 1 16 

7T 

16 



Lt [Lf] 2 + - L 3 , 
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Leading Logarithmic Approximation 




_ol I I I I I I I I I I 

-6 -5 -4 -3 -2 -1 

10 10 10 10 " 10 " 10 10 

w = w* 

Figure 1: Imaginary parts gf_ x of the MHV remainder function in MRK and LLA through 10 
loops, on the line segment with w = w* running from to 1. The functions have been rescaled 
by powers of 4 so that they are all roughly the same size. 



9f\w, w*) = 1 [L - ] 2 [Lif + ^- 6 Lt [Lo ] 4 - ^ L+ [L» ] 2 + 1 [L " ] 2 Ca (4.23) 
"12 [L ^ ]2 Cs + 240 [L ^ ]5 ~ 24 L ° ^ ^ + 384 L * + 8 ^ + 12 ^ 

9?\ W ^) = ~ [L,] 2 [L»f + ^ [L^[Lt] 2 - g [L 2 1 2 + ^ [L 1 4 [^] 2 -^C 3 2 (4-24) 
+^ [^ol 2 M 4 - ^ [^ol 2 [tf] 2 " ^ ] 3 ~§- 4 Lt Li [L -] 2 

"48 ^ [L °" ]2 Cs + 69T^0 [L ° ]6 " 3lo6 [L ° 14 + 4^ [Ls ~ l]2 " ^ L °~ ^ [ ^ ]2 
- ol ^ + ^ 3 " I Ca + ^ [Li? - £ m + ± [Li? + ^ L; L~ 2 

7 1 1 7T 2 __ 9 _ _ 

_ 24 ^ 2 ^ 2 ' 1 ' 1 + 192 ^° ^ 4,:l + 16 ^° ^ 3 ' 2 ~^ 24 ^° ^ 2>1 + 16 ^° 2,1,1,1 

+ 64 5 1 + ^2 1 ^ ~ 48 1 ^ + 32 1 Cs + 12 1 Cs " 32 3 Cs ' 
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Next-to-Leading Logarithmic Approximation 



100 



50 



3^ 



I 

-1 



-50 



-100 



I Ml 1 1 I I I I III 1 1 I I I I II 



/ / v> 




4'-g (2) 

.4 3 -gf 
4 4 gf 

4 6 g< 7) 
4 7 .gf 
4 S -gf 



lO" 3 

w = w* 



_l I I 

-2 -1 

10 10 10 



Figure 2: Imaginary parts gf}_ 2 of the MHV remainder function in MRK and NLLA through 9 
loops. 

sfW) = [L2? [L ] 2 - ^ [Lif [Lj 

29 _ 77 _ ^ 



29 

"9216 ^ 

— [^] 2 
48 L 2,lj 

43 



[l; 



3840 
3 



r]2+ ^ ]2[L " ]4+ ^ [L ; ]4 ^ f 

[L °" ]3 " 5120 ^ Lt [L ° ]2 + 512 Lt [L ° ]2 C3 + 1382400 ^ 

~ ^ L-j [L+] 2 - — L+ [L|] 3 - - [L+] 3 C 3 + ^ [If 

o 29 1 25 3 

960 L * L 2 + 24 L 2 L 2,i,i - ^ L o W,x - ^28 ^ 
301 



(4.25) 



+ 128 " 3 ^48 153 ' 

2?K «>*) = i^SQ ^ [L * ]2 [L « ]2 + i [L2 ~ P [L ^ ]3 " ¥ [L2 " ]2 Cs " 3^ ^ [L °~ ]2 [L ^ ]2 (426) 
+ ^M 2 [^ 2 C3 + ^ra 2 [^ + I ^ 



'2,1,1,1 ^ 
1 

48 



301 L+L+4- 7 L+L+ + 1 L+L+ 3 /V 
3840 Ls Ll + 48 Ll ^ + 12 1 2 < 2 ' x " 128 1 Cs 

^ 2 

S3 ' 



— [Lo ] 4 [L+] 3 - — [L ] 2 [L+ 
18432 576 L J L 1 

"9216 L ^ L * [L °" ]3 " 768 L * [L ° ]2 + 221184 L * [L ° ] ° " 55296 Lt ^ 

"9216 Lt [L » ]A + 2304 L * [L °" ]2 ~ 128 ^ ^ + 768 [L ° ] * Cs 
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-m Lt [Lt]i + is L > ^ ~ rs **» W|2 + k Lt -» [Lt ,2 " s w > 4 Cs 

+ | Wl 2 C= + £ Ml 2 C S + ^ [Li] 7 - ^ [iff 4 I, LI [L t ? 
L 4 L 2 Lf + — — L 2 L Lti — — L 2 L 2 n Lf — — — L L 41 L\ 



192 4 z 1 192 z u a ' L 6 z 1 768 

13 . 5tt 2 + 23 + 21 

_ 384 ' 2 + 144 + 384 ,:L _ 64 

+ 3 L-7- L+ 2157f V 29 L+L+t 19 L+^ 2 + 1 L+ 
+ 16 ^,1,1,1^ -^5I L 5 " 384 ^1 L 3Cs-— ^Cs + ^r 

151 Tj _ 3 ^ 27 ^ 5tt 2 ^ 7 ^ 5tt 2 



-^5,1,1 on 4,1,2 -^4,2,1 , Q "^3,1,1 aA -^3,3,1 70 2 . 2 ,1 



128 32 64 4 ' z,i 48 a ' M 64 ^ 72 

13 + I + 3 + 

+ ^ -^3, 1,1, 1,1 + 2 2,1,2,1,1 + 2 -^2,2,1,1,1 • 

We present only the imaginary parts, as the real parts are determined by eq. (I2.19p . However, 
as a cross-check of our result, we computed the h$ explicitly and checked that eq. (I2.19P is 
satisfied. Furthermore, we checked that in the collinear limit w — > our results agree with the 
all-loop prediction for the six-point MHV amplitude in the double-leading-logarithmic (DLL) 
and next-to-double-leading-logarithmic (NDLL) approximations of ref. 

g^DLLA = in (X (nj -\- HJ* 



1 - J (2^/a\og\w\ log(l - uijj 



h ( 2A/alog \w\ log(l — Mi) ) 

Re (gflNDLLA) = ^ a 3/2^ + w *j | w | V ^ = ^ = / (4.27) 

vAog(i - ui) 

— 7r 2 a 2 (w + w*) log \w\ Iq ^2^/alog \ w\ log(l — uijj , 

where Io(z) and h(z) denote modified Bessel functions. 

Let us conclude this section with an observation: All the results for the six-point remainder 
function that we computed only involve ordinary ( values of depth one for some k), despite 
the fact that multiple ( values are expected to appear starting from weight eight. In addition, 
the LLA results only involve odd £ values - even ( values never appear. 



5 The six-point NMHV amplitude in MRK 

So far we have only discussed the multi-Regge limit of the six-point amplitude in an MHV 
helicity configuration. In this section we extend the discussion to the second independent helicity 
configuration for six points, the NMHV configuration. We will see that the SVHPLs provide the 
natural function space for describing this well. 
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The NMHV case was recently analyzed in the LLA [43]. It was shown that the two-loop 
expression agrees with the limit of the analytic formula for the NMHV amplitude for general 
kinematics (67], and the three- loop result was also obtained. Here we will extend these results 
to 10 loops. 

Due to helicity conservation along the high-energy line, the only difference between the MHV 
and NMHV configurations is a flip in helicity of one of the lower energy external gluons (labeled 

by 4 and 5). Instead of the MHV helicity configuration (+- 1 H — ), we consider (+- 1 1 — )• 

The tree amplitudes for MHV and NMHV become identical in MRK [43]. In this limit, we can 
define the NMHV remainder function -Rnmhv m the same way as in the MHV case (11.11) . 

/i NMHV I _ 4 BDS v ov „/d \ 

^6 | MRK — ^6 X ex Pl-ftNMHVj • \O.L) 

Recall the LLA version^ of eq. f[2~9l : 

At LLA, the effect of changing the impact factor for emitting gluon 4 with positive helicity to 
the one for a negative-helicity emission is simply to perform the replacement 



— n ~> n ( 5 - 3 ) 

-IV + % 

in eq. ( 15. 2p . obtaining [43] 

«f„v - ~ E / ,„J„ (1 - • (5-4) 



2 ^ J-„o (»"+! 

n=—oo v 

The NMHV ratio function is normally defined in terms of as the ratio of NMHV to MHV super- 
amplitudes A, 

„ -4nmhv /_ _x 

^ NMHV = — j • (5-5) 

^•MHV 

However, in MRK, because the tree amplitudes become identical, it suffices to consider the 
ordinary ratio, which in LLA becomes 

4 LLA 

-pLLA _ NMHV ^ rrJ reLLA dLLA \ (r p\ 

' NMHV ~~ a LLA ~~ ^PV-^NMHV rt MHVJ ■ 
MHV 

Thus eq. (15.41) . together with eq. (15. 2p . is sufficient to generate both the remainder function and 
the ratio function in LLA. 



The distinction between R and exp(i?) is irrelevant at LLA, because the LLA has one fewer logarithm than 
the loop order, so the square of an LL term has two fewer logarithms and is NLL. 
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Comparing eq. (15.41) to eq. (12. 9p . we see that in MRK the MHV and NMHV remainder 
functions are related by 

^* ^ LLA 

~R\m\r ■ (<-*-7) 



R 



LLA 
NMHV 



dw 



w dw 



MHV • 



It is convenient to write this equation slightly differently. First, define a sequence of single- valued 
functions f^'(w,w*) in analogy with eq. (E 



R 



LLA 
NMHV 



2m JVlog^l 



Ml j 



1=2 



1 + W* 



f (l) (w,w*) + 



w 



1 + w 



-f {l) (-,- 

* \w w* 



(5.8) 



Then eq. ( I5.7P can be used to relate f^\w,w*) to gf\{w,w* 

d m , 1 



a w u (0 , 
aw — ~ — g\_ x [w, w 



-f 



(0, 



w, w 



w 



\w w 



(5.9) 



w dw* ul 1 v 1 + w 

In Section H] we computed the MHV remainder function in the LLA in the multi-Regge limit up 
to ten loops. Using these results and eq. (I5.9|) . we can immediately obtain NMHV expressions 
through ten loops as well. Indeed, gf^iw^w*) is a sum of SVHPLs, so the differentiation -^L. 
can be performed with the aid of eq. (I3.36p . The result is again a sum of SVHPLs with rational 
coefficients 1/(1 + w*) and w*/(l + w*). As such, the differential equations (I3.36P also uniquely 
determine the result of the w-integral as a sum of SVHPLs, up to an undetermined function 
F(w*). This function can be at most a constant in order to preserve the single- valuedness 
condition. It turns out that to respect the vanishing of the remainder function in the collinear 
limit, F(w*) must actually be zero. 

To see how this works, consider the two loop case. From eq. ( 14. lip . 



(2)/ * 

01 [w,w 



7 [Lt? - 77 [Li 



1 1 1 

2^i,i + 7^0,1 + T^i.o- 



(5.10) 



4 16 

Recalling that (w, w*) = (—z, —z), first use the second eq. (I3.36P to take the w* derivative, which 
clips off the last index in the SVHPL, with a different prefactor depending on whether it is a '0' 
or a '1' (and with corrections due to the y alphabet at higher weights): 



* ^ (2)/ *\ 

W dw^ 91 ( W,W ) 



1 

2 

w 



w ■ 



1 + w* 

1 



Ci- 



w 



4 V 1 



1 + w* 



1+W* 



4 



4 



(5.11) 



4 4 

Next, use the first eq. (13.361) to perform the w-integration. In practice, this amounts to prepending 
a '0' to the weight vector of each SVHPL, 

d 



dw — -^gf ] 
w ow* 



w 



1 + w" 

1 

1 + w* 



■4^0,1-4^0,0 



1 




+ 

1 + w* 


4^0,1 



(2), 



W, W 



+ w* Vw' W*J ' 



(5.12) 



5 Ref. [43] defines a similar set of functions, /;, which are related to ours by ]i = —jf^ 2 \ fa = §/ , etc. 
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where 



w 



4^0,1 

- L 2 + - L L x 



(5.13) 



(log \w\ 2 log(l + w*) - U 2 (-w) + Li 2 (-w*)) 



This result agrees with the one presented in ref. [43]. Furthermore, we can check that the inversion 
property implicit in eq. (15.12)) is satisfied, 



\w w* 



log log 1 + — -Li 2 + Li 2 1 



1 1 



1 



w 



1 



w 



- log 2 \w | 2 — log |w| 2 log(l + w*) + Li 2 (— w) — Li 2 (— w* 



Moving on to three loops, we start with the MHV LLA term, 



g { 2 \w,w*) 



1 



1 



o 3 12 L 1J 

— 777^0,0,1 + 0^0,1,0 + tAu.i + 77^1,0,0 + 7^1.0,1 
16 8 4 16 4 

As before, we can take derivatives and integrate using eq. (I3.36p . 

w* 



w dw 



1 + w* 

1 



+ 



1 + w* 



777 ^0,0,0 — 0^0,0,1 — 777^0,1,0 — 7^0,1,1 
lo 8 lo 4 



1 1 1 

- £0,0,1 + 777 £0,1,0 + 7 £o,l,l 
8 lo 4 



and we find, 

f {3) (w,w* 



77 £0,0,1 + 777 £0,1,0 + 7 £0,1,1 
8 lo 4 



1 

16 



— Lo 1 -|- — Lo ~\~ — Ln Lo -\- Ln Li 

^ z,ii Q 1 u zi Qr) (J 1 

1 



1 

32 



- i log 2 H 2 log(l + w*) + log(-w;)(log 2 (l + w*) - log 2 (l 
+ 2 C2 log |1 + w\ 2 + - log M 2 ( Li 2 (-w) - Li 2 (-to*) 



w 



- 2 log |1 + w\ 2 U 2 (-w) - 2 Li 3 (l +w)-2 Li 3 (l + iw*) + 4 ( 3 



(5.14) 



(5.15) 



(5.16) 



(5.17) 



The last form agrees with the one given in ref. [13] , up to the sign of the second term, which we 
find must be +1 for the function to be single- valued. 



30 



Continuing on to higher loops, we find, 
f^( w , w *) = A Ll ( 3 + l L 2X1 - I L 31 + 1 L \ - 1 L 4 + 1 Lx L 2>1 -^L X L 3 (5.18) 

"4 LojL ' + 4 L ° L2 "Tk LoL3+ 2o6 L ° Ll + lls^^ 

+ 777 Lo -^1 -^2 ? 
lb 

/CB)^, uf) = -1 L 2 C 3 - ^ L C 3 + \ - ~ £2,2,1 + ^ ^4,1 + 7^ ^3,2 (5.19) 

H - — L\ Lo 1 1 ~h Ln Lo 1 1 — Zvi Z/Q 1 -\- Li Lo — La La — Lp. Lo 1 

8 " 16 * ' 16 1 3 ' 32 2 64 1 4 96 ' 

— L? Lo -\- Lt) Ln — Ln La — Ln L? -\- Ln Lo — Lft Li 

96 1 192 2 256 384 1 1152 1536 

~h Lr\ L^ -\- Lfi L\\ — Ln L'i 1 -h — Ln Li Lo i — — Ln Li L^ 

768 1 18432 1 192 3,1 16 ' 48 3 

+ 64 L ° L ' L2 + 768 Ll L2 ~ 8 L ' iX1 ' 

+ 768 Le " 768 ^ 4 ' 2 + 32 ^ + 8 Ll £a « 1 ' 1 ' 1 ~ 16 U U ^ + 16 

-\- Li L '3 9 "I - L, Lo 1 — Lo Li 1 -|- Ln Lo i i — L? Lo 1 -f- L? C.3 

24 ' 32 ' 32 ' 96 " 96 1 ' 96 1 

— L? Lo — Ln Li. i i — Li C's L? Li — Ln Li -\- Li Cn 

128 1 2 192 " 192 s 192 1 256 1 384 s 

— Ln L<\ o — — Ln La i H - Ln 1j\ — Lr\ La ~\- — Ln Lk 

512 3,2 768 ' 960 1 2560 7680 5 
1 . 1 , 5 , 5 _ 5 



— TLn L'i -\- L"r\ Li H- — Lo i C'x ~h L\\ Lx. -\- L?\ Ln 

18432 3 73728 96 ' S 384 5 2048 2 

-\- Ln L? -(- Li La i -)- Ln L? — Ln L '3 i — Lo La 

4096 1 64 ' 1536 1 1536 ' 1536 



11 ~4 T 19 r3~ 1 1 



3 



L n L 2 — L n Lo i -|- — Lo Li Lo i i — — Li L 2 C 

184320 9216 ' 16 " 24 s ' 

1 1 2 1 1 

— L L\ Ls,i + — L L 1 L 2 ,i — — L L,2,i L2 — — L\ L 3 L 2 

h> LlLlL2 ~^2 L ° Ll U + 3^4 L ° LlLl ~^6 Ll Ll Ux 

512 Ll Ll Cs " 96 L ° LlCz ~ 768 U U Cs " 1536 U U U 
11 .0 . , 19 , ,0 . 49 



L Li L3 — — — Lq L\ L3 + Lq Li L 2 . 



2048 u 768 1 18432 

The remaining expressions through 10 loops can be found in computer- readable format in a 
separate file attached to this article. 



31 



6 Single-valued HPLs and Fourier- Mellin transforms 



6.1 The multi-Regge limit in (v,ri) space 

So far we have only used the machinery of SVHPLs in order to obtain compact analytic ex- 
pressions for the six-point MHV amplitude in the LL and NLL approximation. However, this 
was only possible because we knew a priori the BFKL eigenvalues and the impact factor to the 
desired order in perturbation theory. Going beyond NLLA requires higher-order corrections to 
the BFKL eigenvalues and the impact factor which, by the same logic, can be computed if the 
corresponding amplitude is known. In other words, if we are given the functions gb (w,w*) up 
to some loop order, we can use them to extract the corresponding impact factors and BFKL 
eigenvalues by transforming the expression from (w, w*) space back to (is, n) space. The impact 
factors and BFKL eigenvalues obtained in this way can then be used to compute the six-point 
amplitude to any loop order for a given logarithmic accuracy. 

In ref. [2B] the three-loop six point amplitude was computed up to next-to-next-to-leading 
logarithmic accuracy (NNLLA), 

g^(w, w-) = I Li + \ L+ u - i L* \Lt? - | Li \L- f - i Lt L 2 " t L„" 



71 

384 



Lt [Lo? ~ ~ 4 C 5 - y 7" ~ ^Lf? - '-[L^ Lf} 



+ \diC,{[Lif - \\L- f) - \d 2 Li{\Ltf - \[L^} + 1 [Lt 
4 3) (to, «;•) =^LfLt + ± L- %x L - 1 [Lif - 1 [L -] 2 [Lif 



(6.r 



5 



1536 



[L^ + lLfCs 



where d\, d>2 and 7" are some undetermined rational numbers. (To obtain eq. ( 16. ip from ref. [28] 
one also needs the value for another constant, 7' = —9/2, or equivalently 7'" = 0, which was 
obtained in ref. [3D] using the MRK limit at NLLA.) 

These functions can be used to extract the NNLLA correction to the impact factoid. Indeed, 
the NNLL impact factor has already been expressed [3D] as an integral over the complex w plane, 



<m = (-ir ( - 2 + n - ] f —p(w,w*) ki— f - ) 



7T \ W 



2 



6 In principle we should expect the amplitude to NNLLA to depend on both the NNLLA impact factor and 
BFKL eigenvalue. The NNLL BFKL eigenvalue however only enters at four loops, see Section [7^1 
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where the kernel p(w,w*) is related to the three-loop amplitude in MRK, 



p{w,w*) 



(3)/ *\ , i |1 + ^| 2 (3) / 2 \l+w\ 2 

g' Q '{w,w ) + log — — — g\ '(w,w ) + log — = — + vr g\ J (w,w ) 



, , \l+w\ 2 ( . 2 |l + uf 11 
+ log , , C2 log 2 , , - — C4 



\w\ \ \w\ 2 



(6.3) 



However, no analytic expression for $Reg(^' n ) ls y e ^ known. Indeed, an explicit evaluation of the 
integral (16. 2p would require a detailed study of the integrand's branch structure, a task which, if 
feasible in this case, does not seem particularly amenable to generalization. 

Here we propose an alternative to evaluating the integral explicitly. The basic idea is to write 
down an ansatz for the function in (z/, n) space, and then perform the inverse transform to fix 
the unknown coefficients. The inverse transform is easly performed using the methods outlined 
in Section HI so we are left only with the task of writing down a suitable ansatz. To be precise, 
consider the inverse Fourier-Mellin transform defined in eq. (14.41) . Our goal is to find a set of 
linearly independent functions {J^i} defined in (u,n) space such that their transforms {X[J-i]}: 

1. are combinations of HPLs of uniform weight, 

2. are single- valued in the complex w plane, 

3. have a definite parity under Z2 x Z2 transformations in (w,w*) space, 

4. span the whole space of SVHPLs. 

Through weight six, we find empirically that this problem has a unique solution, the construction 
of which we present in the remainder of this section. In particular, we will be led to extend the 
action of the Z 2 x Z 2 symmetry and the notion of uniform transcendentality to (v, n) space. 

6.2 Symmetries in (is, n) space 

Let us start by analyzing the Z 2 x Z 2 symmetry in (u, n) space. It is easy to see from eq. (14.41) 
that 

X[J r (z/, n)](w*, w) = X[J r (z/, —n)](w, w*) , 

X[T(u, n)} (-, -M = T[T(-u, -n)](w, w*) . (6 ' 4) 
\w w* J 

In other words, the Z 2 x Z 2 of conjugation and inversion acts on the (u, n) space via [n -H- —n] and 
[v -H- —v, n -H- —n], respectively. Hence, in order that the functions in (w, w*) space have definite 
parity under conjugation and inversion, J T (u,n) should have definite parity under n <H> — n and 
v — v. Our experience shows that the n- and //-symmetries manifest themselves differently: 
the z/-symmetry appears as an explicit symmetrization or anti-symmetrization, whereas the n- 
symmetry requires the introduction of an overall factor of sgn(n). For example, suppose the target 
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(w Htt>*,WH l/w) 


[y -H- —u, n -H- —n) 


T{y, n) 


(+,+) 


[+,+] 


1/2 \f(v,\n\) + f{- 


-f, 


W)] 


(+,") 


[-,+} 


1/2 [f(u,\n\)-f(- 


-f» 


1*1)] 


(",+) 


[--} 


1/2 sgn(ra) [f(u, |n|) - 


/(■ 


-i/, |n|)] 


(-,") 


[+,-} 


1/2 sgn(n) |n|) + 


/(" 


-i/, |n|)] 



Table 2: Decomposition of functions in (u, n) space into eigenfunctions of the Z2 x Z2 action. Note 
the use of brackets rather than parentheses to denote the parity under (z/, n) transformations. 

function in (w,w*) space is odd under conjugation, and even under inversion. This implies that 
the function in (u, n) space must be odd under n -H- — n and odd under v <H> —v. Such a function 
will decompose as follows, 

F{v,n) = - sgn(n) [f(u, |n|) - f(-v, |n|)] , (6.5) 

for some suitable function /. See Table [2] for the typical decomposition in all four cases. 
Furthermore, in the cases we have studied so far, the constituents f(v, |n|) can always be ex- 
pressed as sums of products of single- variable functions with arguments ±iv + \n\/2, 

f(u, \n\) = J2 c i HfjA$kiv+ \n\/2), (6.6) 

j k 

where Cj are constants, Sf. G {+1,-1}, and the fj,k(z) are single- variable functions that we now 
describe. 

6.3 General construction 

The functional form of J 7 i(u,n) can be further restricted by demanding that the integral (14.41) 
evaluate to a combination of HPLs. To see how, consider closing the //-contour in the lower 
half plane and summing residues at poles with Im(z/) < 0. A necessary condition for the result 
to yield HPLs is that the residues evaluate exclusively to rational functions and generalized 
harmonic numbers, e.g., the Euler-Zagier sums defined in eq. ( 13.101) . This condition will clearly 
be satisfied if the fj,k(z) are purely rational functions of z. Less obviously, it is also satisfied 
by polygamma functions. Indeed, the polygamma functions evaluate to ordinary (depth one) 
harmonic numbers at integer values, 

^(1 +n ) = - lE + Z x [n) and ^(1 + n) = (-l) h+1 k\ (G+i - Z k+1 {n)) , (6.7) 

where ip^ = ip', ip^ = ip", etc. Referring to eq. ( 13.91) . we see that all HPLs through weight three 
can be constructed using ordinary harmonic number^]. 

7 Harmonic numbers of depth greater than one do appear at weight three; however, after applying the stufflc 
algebra relations for Euler-Zagier sums, they all can be rewritten in terms of ordinary harmonic numbers of depth 
one, namely Z 1A {k - 1) = \ Z Y {k - l) 2 - \ Z 2 (k - 1). 
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We therefore expect the fj,k(z) to be rational functions or polygamma functions through 
weight three. Starting at weight four, however, ordinary harmonic numbers are insufficient to 
cover all possible HPLs. Indeed, at weight four, the HPL 



z K 



H h2A (z) = J2jZ 2A (k-l) (6.1 



fc=i 



requires a depth-two suiro, Z 2j ±(k — 1). A meromorphic function that generates Z 2j i(k — 1) was 
presented in ref. [50]. It can be written as a Mellin transform, 



F 4 (A) = M 



Li 2 (x) 



1 — x 

where the Mellin transform M is defined by 



(AO, AeC, (6.9) 



i 

JV 



M[(f(x)) + ](N)= dx(x N -l)f(x). (6.10) 
Jo 

If A is a positive integer, then F^N) evaluates to harmonic numbers of depth two, 

F 4 (A) = Z 2A (N) + Z 3 (N)-( 2 Z 1 (N), NeN. (6.11) 

Going to higher weight, new harmonic sums will be necessary to construct the full space of HPLs, 
and, correspondingly, new meromorphic functions will be necessary to give rise to those sums. 
The analysis of refs. [^9l 150] |5T] uncovers precisely the functions we need_|. They are summarized 
in Appendix [HI Through weight five, three new functions are necessary: F4, F^ a and F-j. 

There is one final special case that deserves attention. Unlike the other SVHPLs, the pure 
logarithmic functions [L^] h diverge as \w\ 0. These functions have special behavior in (z/, n) 
space as well, requiring a Kronecker delta function: 

T\R 1 f lW*M f f + °° dv U„\*» 5n >° t L o] fc+1 (aio\ 

^WHH-LH) {-*) ^Tf H M^ = {kTTy.- (6 ' 12) 

Altogether, we find that the following functions of z = rkiv + \n\/2 are sufficient to construct all 
the remaining SVHPLs through weight five: 

/;.*(*) G |l, p V(l + + z)^"{l + z)rf"{\ + z), Fi(z), F 6a (z), F 7 (z)j . (6.13) 

However, as we will see, not all combinations of elements in the list ( 16 . 1 3[) lead to functions 
of (w,w*) that are both single-valued and of definite transcendental weight. Instead we will 
construct a smaller set of building blocks that do have this property. 



8 Another depth-two sum appears in .#1,1,2(2;) = J^StLi IT^ 1 . 2 ^ — U kut ^ ne * w0 are rented by a stufHe 
identity, Z 2 . 1 (k - 1) + Z 1>2 (k - 1) = Z 2 (k - 1) Z x {k - 1) - Z 3 {k - 1). 

9 Actually, in refs. [3H1 [501 EI] a more general class of functions is denned. It involves generic HPLs that are 
singular at x = — 1 as well as at x = and 1. As wc never encounter these HPLs in our present context, we do 
not discuss these functions any further. 
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6.4 Examples 



Let us see how to use the elements in the list (I6.13P to construct SVHPLs. The simplest case is 
f(v, \n\) = 1. Referring to Tabled only two of the four sectors yield non-zero choices for J- ' . One 
of these, J 7 = sgn(n), produces something proportional to H\ — Hi, which is not single- valued. 
This leaves J 7 = 1, which should produce a function in the (+, +) sector. Closing the //-contour 
in the lower half plane, and summing up the residues at v — — i\n\/2, we obtain the integral of 
eq. (USD, 



J[l] = 2L\ 



(6.14) 



indeed a function in the (+, +) sector. Including the special case Lq from eq. (I6.12p . this 
completes the analysis at weight one. 

The next simplest element is 1/z, yielding f(v, |n|) = l/(iv + \n\/2). It generates two single- 
valued functions, one in the (+, — ) sector and one in the (— , — ) sector (using the (w, w*) labeling 
in the first column of Table |2]). Symmetrizing as indicated in Tabled the two functions in (v, n) 
space are T = — V and J 7 = N/2, with the useful shorthands 



V 



1 



N = sgn(n) 



1 



i l n l 



i \ n \ 



iv 



-iv + y 



|n| 2 
4 

n 



(6.15) 



M 2 



The transforms of these functions yield two of the four SVHPLs of weight two. 



X[V] = -L " Lf 
I[N] =4L 2 ". 



(6.16) 



A third weight-two function is the pure logarithmic function [Lq ] 2 , a special case already con- 
sidered. To find the fourth weight-two function, we turn to the next element in the list (I6.13p . 
^(1 + z). On its own, it does not generate any single- valued functions; however, a partic- 
ular linear combination of {1, 1/z, ^(1 + z)} indeed produces such a function. Specifically, 
f(v, \n\) = 2*0(1 +iv+ \n\/2) + 2-y E - l/(iv + \n\/2) generates the last weight-two SVHPL, 
which transforms in the (+, +) sector. The function in (v, n) space is actually the leading-order 
BFKL eigenvalue, E vn , 



\n\ 



1+W+y | + ^ 



\n\\ sgn(n)N 



,n ■ 



(6.17) 



and its transform is the last SVHPL of weight two, 



I[E v>n ] = [LJY - - [Lq] a ■ 



(6.11 
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The next element in the list (I6.13P is + z). Like + z), + z) does not by itself 
generate any single- valued functions; however, there is a particular linear combination that does, 
and it is given by f(v, \n\) = 2ip'(l + iv + |n|/2) + 1/ (iv + \n\/2) 2 . Notice that, for the first time, 
the product in eq. (16. 6p extends over more than one term (in this case, = j\ : 2 = 1/ (if + M/2), 
but in general the fj t k will be different). The function in (v, n) space is, 

T = iff' ( 1 + iv + H ) - V' (l-iv + ^f)- sgn(n)NV = D v E v , n , (6.19) 



2 / V 2 , 

where D v = —id u = — z djdv. The main observation is that the basis in eq. ( 16 . 1 3[) can be modified 
to consistently generate single- valued functions: 1/z is replaced by V and TV, -0 is replaced by 
E VtTl , and ip^ is replaced by D^E v ^ n . 

Furthermore, as mentioned previously, the basis at weight four requires a new function F 4 (z) 
that is outside the class of polygamma functions. Like the polygamma functions, F 4 (z) does 
not by itself generate a single-valued function; it too requires additional terms. We denote the 
resulting basis element by F 4 . It is related to the function F 4 (z) in eq. (16.91) by, 



~ / \fl\ \ { \Tl\ \ 

F 4 = sgn(n) I F 4 (iv + U J + F 4 ( - iv + U J - -D^, n - -iV 2 ^„ - - 



(6.20) 



where 



V>- = ^(l + ^ + y) ~ H 1 ~ zV + t) ' (6 ' 21) 

Appendix [B] contains further details about the functions in (v, n) space, including the basis 
through weight five and expressions for the building blocks F Sa and F 7 generated by the functions 
F 6a (z) and F 7 (z). 

Finally, we describe a heuristic method for assembling the basis in (v, n) space. The idea is 
to start with the building blocks, 

{l,N,V,E v>n ,F 4 ,F 6a ,F 7 }, (6.22) 

and piece them together with multiplication and ^-differentiation. These two operations do not 
always produce independent functions. For example, 

D U N = 2NV and D V V = ^N 2 + V 2 . (6.23) 

The building blocks have definite parity under v •H- — v and n -H- — n which helps determine 
which combinations appear in which sector. Additionally, we observe that they can be assigned 
a transcendental weight, which further assists in the classification. The weight in (w,w*) space 
is found by calculating the total weight of the constituent building blocks in (v, n) space, and 
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weight 


(v —v, n -H- —n) 


weight 


[y -H- — v, n <H- — n) 


1 





[+,+] 




1 


[+,+] 


D v 


1 


K +] 




3 


[+,-] 


V 


1 


h,+] 




4 


h,-] 


N 


1 


[+,-] 


F 7 


4 


K +] 



Table 3: Properties of the building blocks for the basis in (y, n) space. 



then adding one (to account for the increase in weight due to the integral transform itself). The 
relevant properties of the basic building blocks are summarized in Table [3j 

As an example, let us consider the function ND v E v ^ n . Referring to Tabled the transcendental 
weight is 1 + 1 + 1 — 3 in (u, n) space, or 3 + 1 = 4 in (w, w*) space. Under [u -H- —u, n -h- —n], 
N has parity [+,—], D v has parity [—,+], and E u>n has parity [+,+], so ND u E v ^ n has parity 
[—,—]. We therefore expect this function to transform into a weight four function of (w,w*), 
with parity (— ,+) under (w -H- w*,w -H- 1/w) (see Table |2]). Indeed this turns out to be the 
case. A complete basis through weight four is presented in Table HJ 



7 Applications in (z/, n) space: the BFKL eigenvalues and 
impact factor 

7.1 The impact factor at NNLLA 

In this section we report results for g^' and an d discuss how to transform these functions to 
(z/, n) space using the basis constructed in the previous section. We then give our results for the 

(2) (3) (9) 

new data for the MRK logarithmic expansion: $j^ g , $R,eg' an< ^ ^,n- 

Before discussing the case of the higher-order corrections to the BFKL eigenvalue and the 
impact factor, let us review how the known results for E VtTl , El]l and $j^ g fit into the framework 
for (u,n) space that we have developed in the previous section. First, we have already seen in 
Section [6] that the LL BFKL eigenvalue is one of our basis elements of weight one in (u, n) space 
(See Table |3]). Next, we know that the first time the NLL impact factor $Reg appears is in the 

(2) 

NLLA of the two-loop amplitude, g$ (w,w*), which is a pure single- valued function of weight 
three. Following our analysis from the previous section, it should then be possible to express 
^Rog as a P ure function of weight two in (r/, n) space with the correct symmetries. Indeed, we 
can easily recast eq. ( I2.17P in terms of the basis elements shown in Table [31 



0^) = -^,n-^ 2 -C2. (7.1) 

Similarly, the NLL BFKL eigenvalue can be written as a linear combination of weight three of 
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Z9 X Z9 


(7/; 7/;*l basis 


(7/ r) ] basis 


dimension 


1 


( 1 1 \ 

( + , +J 


^1 


1 
1 


1 
1 


( ~\~ ) 


r — 
L 




1 
1 


f i\ 

\ ' ' ) 






n 

u 











2 


( _i_ 1 ^ 

1 + 7 + ) 


r r +12 rr-12 
1^1 J > i L oJ 






{ | ) 




V 


1 


V 1 'J 






n 




^2~ 


N 


1 


3 


( _i_ 1 *\ 
\ 1 ) + J 


r r +13 [r-l2r+ r + 


T/2 Ar2 7712 


q 



z' 1 ^ 

i+> _ j 


r 7 — 13 r-[r+l2 r- 
L-^o J > ^0 L-^l J > "^2,1 




Q 
O 




r — r — 
-^0 -^2 


T/ AT 


1 
1 








1 


4 


(+,+) 




<WH 3 > £ 3 n , iV 2 ^, n , 
^ /2 -E i /,n, VD v E^ n , D 2 V E V ^ 
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La [L x ] 3 , [L ] 3 L x , L L 3 , L x L 21 


V 3 , N 2 V, VE 2 >n , E„ >n D v E v , n 
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(->+) 


La Li L 2 , I/ 3jl 
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Table 4: Basis of SVHPLs in (w,w*) and (1/, n) space through weight four. Note that at each 
weight we can also add the product of zeta values with lower-weight entries. 



the basis elements in Table [3j 



"4 D u E u<n + - V D V E V)U — ( 2 E Utn — 3 C3 • 



(7.2) 



This completes the data for the MRK logarithmic expansion that can be extracted through two 
loops. 

Now we proceed to three loops. By expanding eq. (14.11) to order a 3 , we obtain the following 
relation for the NNLLA correction to the impact factor, $^ g (z/, n), 

X 

- A<K*gf\w, w*) Lf + ^Lf {-45 [L^f + 120 [L+] 2 + 22 vr 2 } . 



2 (7-3) 

2,(2), 



This expression is exactly 2 p(w,w*), where p was given in eq. (16. 3p and in ref. [ID]. (The factor 
of two just has to do with our normalization of the Fourier- Mellin transform.) 
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To invert eq. (17. 3p and obtain & R ^ g (v,n), we begin by observing that the right-hand side is 
a pure function of weight five in (w, w*) space. Moreover, it is an eigenfunction with eigenvalue 
(+, +) under the Z 2 x Z 2 symmetry. Following the analysis of Section [61 and using the results 
at the end of Appendix [BJ we are led to make the following ansatz, 



,(2) 
Rcg\ 



ai E% n + a 2 N 2 El >n + a 3 + a 4 V 2 E\ n + a 5 N 2 V 2 + a 6 V* 

[D V E V ^ + a 9 E v , n DlE v , n + a w F 4 N (7.4) 
+ an C,2E 2 vn + «i 2 C2^V 2 + a 13 ( 2 V 2 + a u CsK,n + a 15 (3 [S n ,o/ + «ie C4 • 

The aij are rational numbers that can be determined by computing the integral transform to 
(w, w*) space of eq. (I7.4p (see Appendix [B]) and then matching the result to the right-hand side 
of eq. (Q. We find 



^Rcg 
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4 360 



5vr 2 ^ 2 
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C3 E V i 



N 4 



7T 
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(7.5) 



+ 7" 



7T 



6 



£ 2 



Here di, g?2 and 7" are the (not yet determined) rational numbers that appear in eq. (16. ip . We em- 
phasize that the expression for n ) does n °t involve the basis element N F 4 (see eq. ( IB. 52ft ). 

(2) 

That is, v,n) can be written purely in terms of ip functions (and their derivatives). 

To determine the six-point remainder function in MRK to all loop orders in the NNLL ap- 
proximation, we must apply some additional information beyond ^^ g (u,n). In particular, at 

(2) 

four loops and higher, the second-order correction to the BFKL eigenvalue, Ev,n, is necessary. 

In the next section, we will show how to use information from the symbol of the four-loop re- 

(2) 

mainder function to determine ElA- We will also derive the next correction to the impact factor, 
$Reg( v ' n )) w hich enters the N 3 LL approximation. 



7.2 The four-loop remainder function in the multi-Regge limit 

In order to compute the next term in the perturbative expansion of the BFKL eigenvalue and the 
impact factor, we need the analytic expressions for the four-loop six-point remainder function in 
the multi-Regge limit. In an independent work, the symbol of the four- loop six-point remainder 
function has been heavily constrained [52]. In ref. [52] the symbol of is written in the form 

113 

S(4 4) ) = 5>^, (7.6) 

i=l 

where a.i are undetermined rational numbers. The Si denote integrable tensors of weight eight 
satisfying the first- and final-entry conditions mentioned in the introduction, such that: 
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1. All entries in the symbol are drawn from the set {itj, 1 — itj, yi}i=x,2,3, where the y^s are 
defined in eq. (11.41) . 

2. The symbol is integrable. 

3. The tensor is totally symmetric in Ux, u 2 , u 3 . Note that under a permutation Ui — > u a ^, 
a G S 3 , the yi variables transform as yi — > l/y^. 

4. The tensor is invariant under the transformation yi — > 1/yi. 

5. The tensor vanishes in all simple collinear limits. 

6. The tensor is in agreement with the prediction coming from the collinear OPE of ref. [25] , 
We implement this condition on the leading singularity exactly as was done at three 
loops [28] . 

In Section HI we presented analytic expressions for the four-loop remainder function in the LLA 
and NLLA of MRK. We can use these results to obtain further constraints on the free coefficients 
«i appearing in eq. (17. 6p . In order to achieve this, we first have to understand how to write the 
symbol (17. 6p in MRK. In the following we give very brief account of this procedure. 

To begin, recall that the remainder function is non-zero in MRK only after performing the 
analytic continuation (12. 6p . u\ — > e~ 2nl \ui\. The function can then be expanded as in eq. (12. 7p . 

3 

4 4) |mrk = 2th ^"(i [g [ n ) (w J w*)+2mh < n i \w,w*)} . (7.7) 

n=0 

The symbols of the imaginary and real parts can be extracted by taking single and double 
discontinuities, 

3 



n=0 

113 

= -2ni ^ a i A «i ( S i) Imrk 

i=l 

(2mf J2 S H n (l ~ ui) h^(w,w*)] = 5(A^4 4) )| MRK 



(7- 



n=0 
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-2ixif A^ (Si) |mrk 



i=i 



where the discontinuity operator A acts on symbols via, 



a 2 ® . . . <g> a n , if a x = Ux , 

A Ul (a!®a 2 ® ...®a n ) = { (7.9) 

, otherwise. 

5 (o 3 ® . . . ® On) , if ax = a 2 = Ux , . . 

A; (fl!®a 2 ® ...®a n ) = < n . (7.10) 

(J , otherwise. 
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As indicated in eq. (17. 8p . we need to evaluate the symbols Si in MRK, which we do by taking 
the multi-Regge limit of each entry of the symbol. This can be achieved by replacing w 2 and u% 
by the variables x and y, defined in eq. (12. 3p (which we then write in terms of w and w* using 
eq. (I2.4p ). while the y^s are replaced by their limits in MRK [28] . 

1 + w* w*(l + w) 

1 + w w(l + w*) 

Finally, we drop all terms in A^ (Si), k = 1,2, that have an entry corresponding to u\, y±, 1 — u<i 
or 1 — U3, since these quantities approach unity in MRK. In the end, the resulting tensors have 
entries drawn from the set {1 — U\, w, w*, 1 + w, 1 + w*}. The 1 — Ui entries come from factors 
of log(l — Ui) and can be shuffled out, so that we can write eq. (17. 8p as, 

3 113 7 



Y,S[log n (l -u x )]mS[g^(w,w*))] = £ £a,S [log"(l - Ul )\ uiG v „ 

n=0 i=l n=0 

3 113 6 

J2S[log n (l-u l )}mS[h^(w,w*))} = ^5>*<S[log"(l- Wl )]in# i)n , 



(7.12) 



n=0 i=l n=0 



for some suitable tensors Gi tU of weight (7 — n) and H in of weight (6 — n). The sums on the 
right-hand side of eq. (I7.12p turn out to extend past n = 3. Because the sums on the left-hand 
side do not, we immediately obtain homogeneous constraints on the «j for the cases n = 4,5, 6, 7. 
Furthermore, since the quantities on the left-hand side of eq. (I7.12p are known for n = 3 and 
n = 2, we can use this information to further constrain the Oj. Finally, there is a consistency 
condition which relates the real and imaginary parts, 

hf(w, w*) = g*i\w, w*) + ^ gf\w, w*) L\ -\ gf(w, w*) L\ - gf\w, «,*) gj?\w, w*) , 
hf(w, w*) = \ g{ 4) (w, w*) + rr 2 gf\w, w*) - rr 2 gf\w, w*) Lf - ± gj?\w, w*) Lf 

+ y *fV ™*) [Li? + y 2 ^ «»*) ^ + g [^o ] 2 [tf] a " ^ [^o ] 4 

+ JL^ [ L -]2 _ 1^2 [L +14 _ _3_ 4 [L+] 2 _ I [ (2) ( * )]2 

640 1 J 96 1 1 J 160 1 1 J 2 [y ° v ' ;J 

(7.13) 

In total, these constraints allow us to fix all of the coefficients ati that survive in the multi-Regge 
limit, except for a single parameter which we will refer to as a . 

The results of the above analysis are expressions for the symbols of the functions g^ and g^ . 
We would like to use this information to calculate new terms in the perturbative expansions of the 
BFKL eigenvalue u(u,n) and the MHV impact facator $ Reg (^, n). For this purpose, we actually 
need the functions g[^ and g^\ and not just their symbols. Thankfully, using our knowledge of 
the space of SVHPLs, it is easy to integrate these symbols. We can constrain the beyond-the- 
symbol ambiguities by demanding that the function vanish in the collinear limit (w,w*) — > 0, 
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and that it be invariant under conjugation and inversion of the w variables. Putting everything 
together, we find the following expressions for and 

9?(w,w*) = JL [L2 ] 2 [L ~] 2 - | [L 2 ] 2 [Ltf + ^ [L -] 2 [L+] 4 + ^ W W 
+ gg [Lo ] 3 " 771 ^ L 3 M 2 - - g [^i] 2 ~^Lj [Lff 



[^ 2 1 2 [Lo? + ^ [^o] 2 + ^ [^o] 4 [^] 2 + 
1581 j.,^,, „ 9823 _ ,« 871 , , r ^ l0 709 



^ Kl 2 + [iol 6 - ^ ^ i„ [Lt? - ^ Li [Ltf (7-14) 



16 1 ° L UJ 1152 4 
9993 

+ L+ L+ - 157 [L 2 ] 2 [L+] 2 - 256 [L 2 J 2 + 1593 [L+f ( 3 

+ 681 [L+] 2 - 1606 L4 Lo; + 512 L 2 L 2ll - 3371 Ll x 

- 1730 L L7 ;2 - 299 L + 2127 L+ L+ M + 744 Lf L+ 21 

+ 5489 L+ Cs + 256 C 3 2 J + «i tt 2 #f ) (iy , w*) + a 2 n 2 (w, w*) 

+ a 3 7r 2 [gj 2 ^ (w, w*)] 2 + a 4 7r 2 fa 2 (w, w*) + a 5 7r 2 /iq 3 - 1 (w, w*) 
+ a 6 vr 4 gf ] (w, w*) + a 7 ( 3 g£> (w, w*) + a s ( 3 g { 2 ] (w, w*) . 



(4)/ * 

g y '{w,w 



I Lt [L^] 2 [L - ] 2 - i [L 2 ] 2 [L+] 3 - | L+ [L -] 2 [L+] 2 + ^ [L ~] 4 [L+] 3 



+ Jo [L °" ]2 ^ " 192 



^2,1 ^1" [-^0 ] 3 



129 
~64~ 



2^ ^3 [-^0 ] 4 + 32 ^•V^ 1 ^° ] 2 



1 

16 



Lt [Lo? 
7 



1007 
46080 



Li [L~ 



-16 



^£2,1 [-^0 ] 2 + T7 [-^0 ] 2 Cs 



16 

1 L - L 2 - x [L+] 3 + ^ L+ [L+] 2 - 1 Li [Lif + ? - L+ X)1 [L+] 2 



16 

+ 12 [ ^ ]4C3 + 210 



16 

\Li\ 7 



48 

5 1 
L 4 L 2 L 1 — — L 2 L L 31 + -L 2 L 2 l l L 1 



1 - - + 1 - - + 
~*~ 4 ^° "^ 4 > 1 ^ 1 _ 8 ^° ^' 2 ' 1 ^ 3 

l^_i_ 11 3 ^ 

+ 2 4 ' 1 ' 2 "4" A ,2,\ + 7 -^3,3,1 



16 

^ -^0 -^2,1,1,1 ^l" + 2 -^l" ^3 ~ 

25 



125 



±r+ _ _ r 

r, "^2,1,2,1,1 r, 



2,2,1,1,1 



^ Ct + 5 



(7.15) 
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+ — # [Lt? - — Lt [Lff + 205 [L~ 2 f [Lff + 1911 Lf [L^f 



+ 2130 L+ 2>1 [L ~] 2 - 2623 [L ~] 2 C 5 + 992 Lf [L^f + 63 L+ M [Lt] 2 

- 288 Lj 21 [L+] 2 + 2396 [L+] 4 Cs + 1830 Lf [Lf] 2 - 1612 L± L 2 Lf 
+ 1344 L 2 L Lf A - 520 L 2 L 21l Lf + 11839 Lq L^ Lf 

+ 4330 Lq Z£ 2 Lf + 3780 L L 21 Lf + 562 Lq L 2lll Lf 

+ 3556 Lf d + 2256 Lf - 164778 Lj M - 33216 L+ 12 - 89088 L+ 21 

- 33912 L+ 31 - 12048 L+ 2 2 - 17820 - 2928 Lf 121l 



+ 63 C2 g? (w, w*) (w, w*) + W C2 9V (w, w*) g[ z> (w,w*) + b 5 ( 2 h™ (w,w*) 
+ 6 6 C 2 hf (w, w*) + b 7 C 2 g [ o ] (w, w*) + b 8 C 2 (w, w*) + b 9 ( 2 gf (w, w*) 
+ 6 10 Cs hf (w, w*) + b n Cs ^ K w*) + 6 12 Cs K Of 

+ 613 C3 #3 O, w*) + 614 C3 g? (w, w*) + 615 C4 g 2 ] (w, w*) + 6ie C4 fif } O, W*) 
+ 617 C3 C2 g? (w, W*) + 6 18 Cs #1 (u>, U>*) • 



In these expressions, aj for i — 0, . . . , 8, and bj for j = 1, . . . , 18, denote undetermined rational 
numbers. The one symbol- level parameter, a , enters both gf and gf. We observe that enters 
these formulae in a complicated way, and that there is no nonzero value of that simplifies the 
associated large rational numbers. We therefore suspect that ao = 0, although we currently have 
no proof. The remaining parameters account for beyond-the-symbol ambiguities. We will see in 
the next section that one of these parameters, 61, is not independent of the others. 

7.3 Analytic results for the NNLL correction to the BFKL eigenvalue 
and the N 3 LL correction to the impact factor 

Having at our disposal analytic expressions for the four-loop remainder function at NNLLA and 
N 3 LLA, we use these results to extract the BFKL eigenvalue and the impact factors to the same 
accuracy in perturbation theory. We proceed as in Section [7TT| i.e., we use our knowledge of the 
space of SVHPLs and the corresponding functions in (z/, n) space to find a function whose inverse 



- 8784 L 



2,2,1,1,1 



23796 (A+h C2 [L 2 ] 2 Lf + b 2 (2 [L^] 2 Lfgf\ W) w*) 
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Fourier-Mellin transform reproduces the four-loop results we have derived. 

Let us start with the computation of the BFKL eigenvalue at NNLLA. Expanding eq. fl4.ip 
to order a 4 , we can extract the following relation, 



Z[El 2 l] =12 {lLf] 2 + n 2 }gi 4 \w,wn-8Ltg { 2 4) (w,w*)+4g { *\w,w*) 
- 8 Lf tt 2 g® {w, w*) + 2 vr 2 gf\w, w*) [L+] 2 



(7.16) 



1 



-1 



The right-hand side of eq. (I7.16P is completely known, up to some rational numbers mostly 
parametrizing our ignorance of beyond-the-symbol terms in the three- and four-loop coefficient 
functions at NNLLA. It can be written exclusively in terms of SVHPLs of weight six with 
eigenvalue (+, +) under Z2 x Z2 transformations. The results of Section [6] then allow us to write 
down an ansatz for the NNLLA correction to the BFKL eigenvalue, similar to the ansatz ( I7.4p 
we made for the NNLLA correction to the impact factor, but at higher weight. More precisely, 
we assume that we can write E^l = ^ i a i P i , where a« denote rational numbers and Pi runs 
through all possible monomials of weight five with the correct symmetry properties that we can 
construct out of the building blocks given in eq. (I6.22p . i.e., 



P t e El C 2 V D u E v>n , E„ n N F 4 , C 5 , 



(7.17) 



The rational coefficients a.i can then be fixed by inserting our ansatz into eq. (I7.16P and performing 
the inverse Fourier-Mellin transform to (w, w*) space. We find that there is a unique solution for 
the ai, and the result for the NNLLA correction to the BFKL eigenvalue then takes the form, 

42 = $ Ro>> n) - E v>n sg^i/, n) + | D 2 u E u>n £ 2 n + A tf» D 2 V E^ + l -V 2 D 2 v E^ n 



48 v u ' n 



--VDZE Vtn 
- 2 1 



7T 2 3 2 

^ D v E v,n - ^ D v E VjTl V E V1 



7T 



— D v E y ^ n V + - E v<n [D v E Vin ] 

llTT 



^E Utn N 2 + ^E u , n V 2 - 
16 4 



16 
1 

" 8 
35 

16 



D u E u , n N 2 V 
5tt 2 

EL + —EI 



6 



N 2 Cz 



V 2 Cs 



+ —7— C3 + 10 C5 + a S 5 + ^ a i C2 £3,1 + a & U £2 + ^2 a { ( 3 £ lfi 



i=l 



i=7 



(7.18) 



where the quantities S^i, S2, and £n capture the beyond-the-symbol ambiguities in g± \ and £< 
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corresponds to the one symbol-level ambiguity. They are given by, 



+ N * £ 3„ _ ^ V« B „»„ _ | ^ ^^i^^C^i^^ C3 (7 , 9) 

"1 2 771 77i2 I 7~1 Z7 1 Z7>2 vcn n Z7 1 T/3 i O/IO 7Ti [n P 12 



62 DX,n + 70 D v E u , n V E^ n - 760 D u E v>n V + 248 E v , n [D v E v , n ] 

5 • 



-7431 _ Cs - 97 n A^ 2 V 2 + 16072 £ 



£3,1 = ~E u , n N 2 -DlE U:n + 5El n + QE^ n V 2 -2E^ n TT 2 + 8( 3 , (7.20) 
£3,2 = £j„, (7.21) 

£3,3 = I E Ujn N 2 - 3 V + 3£j n + 12 C 3 , (7.22) 

1 q 3 3 2^ 

£3,4 = -~DlE v>n + ~D v E v , n V-~E v , n N 2 --E v , n V 2 --(3-2El n , (7.23) 

£3,5 = \E Vtn N 2 - Z -El tn , (7.24) 
£ 2 = 90£„, n , (7.25) 
5 1)7 = ^ 2 n -^iV 2 , (7.26) 

£1,8 = ^,„. (7-27) 

We observe that the most complicated piece is £5. It would be absent if our conjecture that 
a = is correct. Some further comments are in order about eq. (17. 18ft : 

1. In ref. (10] it was argued, based on earlier work [SHI EH EDI Ej, that the BFKL eigenvalue 
should vanish as (u, n) — > to all orders in perturbation theory, i.e., w(0, 0) = 0. While 
this statement depends on how one approaches the limit, the most natural way seems to 
set the discrete variable n to before taking the limit v — > 0. Indeed in this limit E VtTl 
and Eu)i vanish. However, we find that E^) n does not vanish in this limit, but rather it 
approaches a constant, 

lim<o = -^ 2 C 3 . (7.28) 

We stress that the limit is independent of any of the undetermined constants that parame- 
terize the beyond-the-symbol terms in the three- and four-loop coefficients. While we have 
confidence in our result for E v 2 ) n given our assumptions (such as the vanishing of g£) and 
hn as w — > 0), we have so far no explanation for this observation. 

2. While the (v, n)-space basis constructed in Section [6] involves the new functions F4, Fq g and 
F 7 , we find that Ev,n is free of these functions and can be expressed entirely in terms of ijj 
functions and rational functions of v and n. Moreover, the ip functions arise only in the 
form of the LLA BFKL eigenvalue and its derivative with respect to v. We are therefore 
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led to conjecture that, to all loop orders, the BFKL eigenvalue and the impact factor can 
be expressed as linear combinations of uniform weight of monomials that are even in both 
v and n and constructed exclusively out of multiple ( values^ and the quantities N, V, 
E v<n and D v defined in Section [61 

We now move on and and extract the impact factor at N 3 LLA from the four-loop amplitude 
at the same logarithmic accuracy. Equation (14. ip at order a 4 yields the following relation for the 
impact factor at N 3 LLA, 



X 



(4), 



ALtg?\w,w*)+Ag^(w,w*) + S^g^(w,w*) 

r 2 



,( 4 )/ 



.2 „(3), 



+12 



- 4 L+tt 2 ^ 3) (w,w* 



7T 



(2)/ * 



(7.29) 



In order to determine $^ (v, n), we proceed in the same way as we did for Eb%, i.e., we write 

down an ansatz for $Rg g (^, n) that has the correct transcendentality and symmetry properties 
and fix the free coefficients by requiring the inverse Fourier-Mellin transform of the ansatz to 
match the right-hand side of eq. (I7.29p . Building upon our conjecture that the impact factor can 
be expressed purely in terms of ip functions and rational functions of v and n, we construct a 
restricted ansat z0 that is a linear combination just of monomials of ( values and N, V, D v and 

(2) 

E u n . Just like in the case of E„ t n, we find that there is a unique solution for the coefficients in 
the ansatz, thus giving further support to our conjecture. Furthermore, we are forced along the 
way to fix one of the beyond-the-symbol parameters appearing in g^\ 

15 3 3 9 1 , 1 r. 3 j 1. 1 , 3 j 3 _ on . 

b\ = a\ a 2 a4 H as H 03 H — o 4 o 5 o« H — o 7 H o« H . (7.30) 

8 16 32 16 64 8 16 32 4 32 16 V ' 

The final result for the impact factor at N 3 LLA then takes the form, 



7T 4 7T 4 73-7T 6 

w, w*) \Lf} 2 + —Lf [Ln ] 2 - — [L+] 3 - — 
' ;LlJ 8 lLoJ 3 LlJ 1260 



Lt-2Lttt 
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<^ n) = \ n)] - Efl E v>n - n) £ 2 n - - [DlE v>n 
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(7.31; 
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- D u E v ^ n V D u E u , n — — D u E Utn D v E U:n + - D u E Utn E un — — Ey <n N D v E UjT 
4 24 o o2 



37tt 2 
96 

llTT 2 



E u n D E v n — — D E v n (3 — 7 D v E v n V E v + — 
24 4 16 



£" iV 2 V 



f D u E un E u>n V + 7 D u E un V (3 + 4 [A^,n] 2 £ 2 n - |r iV 2 [£,£„,n] 2 
24 4 16 64 



10 Note that we can not exclude the appearance of multiple £ values at higher weights, as multiple £ values are 
reducible to ordinary £ values until weight eight. 

11 We have constructed the full basis of functions in [y, n) space through weight six and the explicit map to 
(w, w*) functions of weight seven. It is therefore not necessary for us to restrict our ansatz in this way. It is, 
however, sufficient, and computationally simpler to do so. 
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- i V 2 [D„.E„„] 2 + ^ [C„B„,„] 2 + ^ JV a £ 2 „ + A JV 2 V 2 £ 2 „ 

-^2g i^,n 22 u,n 1 ^2 v,n 1 g u,n 

+ l Evn N 2 { 3 -—E„ n ( 3 -—N e -—N 4 V 2 + —N*--N 2 V 4 
16 ' 2 ' s 1024 128 512 16 
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128 24 8 2 S3 11340 ' S ' S 

15 5 8 

+ — iV 2 E 4 n + a V 6 + ^ai ( 2 V a ,i,i + a 6 Q V a ,2 + 2^ ai ^ 3 

i=l i=7 
9 14 16 

+ J] bit 2 v bAi + ^ &i c 3 n, 3> i + 5^ & * c 4 n,2,i + 617 C2C3 n,i,i + h 8 c 5 n,i, 2 , 

i=2 i=10 i=15 

where 7\j r .. parametrize the beyond-the-symbol terms in the four-loop coefficient functions, and 
V% parameterizes the one symbol-level ambiguity, 

-~ E v , n D 4 v E v , n - 1^ DX„ Ca + ^ iV 2 [^,n] 2 + ^ F 2 [A,^] 2 

+^ ^ El n + ^N 2 V 2 £ 2 , n + ^ TV 2 ^ + # >f| + ^ iV 2 C 3 

+ 13777 ^ y2 ^ + 111 ^ + 345 ^ y2 _ 3g5 y _ 3Q 

+ 16 D 2 v E U)n E% n - 420 D u E Utn V E% n + 7 D v E V)n E v>n N 2 V — 760 £ v>n ^ 3 

»5 



-22606 D v E Vtn V (3 ~ 34 [D,^] 2 -E 2 ^ + 1140 V 4 E 2 n + 15231 £ 3 n ( 3 + 6548 E v>n ( 5 



+46992 C 3 2 , 

^a,4,l = o -Ei/.n D 2 E U)n — - D v E v<n E VjU V + — [-D^-E^] 2 — — iV 2 i? 2 n (7.33) 
o 2 o o<i 

-™V 2 E 2 n --Et n + — N 4 + —N 2 V 2 + — N 2 + n 2 E 2 -68E vn ( 3 , 

q o V ri 

V aA2 = --E u , l DlE^ n + -D u E^ n E^ n V + -[D u E^ n ] 2 --N 2 El n (7.34) 
16 4 16 64 



.^y 2 E 2 n --^ 4 n -— iV 4 + -iV 2 ^ 2 -12^ n C3, 
16 v ' n 4 256 64 



^,4,3 = -^.n^X^P^^ (7.35) 



+^ TV 2 V^ 2 + 6 D„E v , n E v , n V - 48 Cs , 
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n,4,9 = 24 (7-48) 

1 3 11 1 25 

n, 3 ,io = --£>Xn + ^D v E u , n V - -El n - -E^ n N 2 - -E„, n V 2 - -( 3 , (7.49) 

n, 3 ,n = ^^iV 2 -!^, (7.50) 

n, 3 ,i2 = -i J D^,„^ + ^J n + i^ in iV 2 + 2C3, (7.51) 

n, 3 ,i3 = i<„, (7.52) 

n,3.14 = -^Xn + ^'n-i^.n^-^^+^Cs + ^.n^ 2 , (7.53) 

o 8 3 3 

n,2,i 5 = (7.54) 

n,2,i6 = ^, n -iiV 2 , (7.55) 

= ^,n, (7.56) 

n,i,2 = (7.57) 



Again, the undetermined function at symbol level, Vq, is the most complicated term, but it would 
be absent if a = 0. 

Finally, we remark that the v — > behavior of $^ g (^, n) is nonvanishing, and even singular 

for £ = 2 and 3. Taking the limit after setting n = 0, as in the case of Ei,,n, we find that the 
constant term is given in terms of the cusp anomalous dimension, 

Jm^M) ~ it + Q^), (7.58) 

2 (3) 

Jim^tM) ~ ^ + ^ + 0(S), (7.59) 

4 (4) 

hm<M) ~ -^2+^ + ^)- (7-60) 

This fact is presumably related to the appearance of 7a- (a) in the factors u a b and S, which carry 
logarithmic dependence on as it) — > 0. It may play a role in understanding the failure of E„q 
to vanish as v — > in eq. f)7.28p . 

8 Conclusions and Outlook 

In this article we exposed the structure of the multi-Regge limit of six-gluon scattering in planar 
M = 4 super- Yang-Mills theory in terms of the single- valued harmonic polylogarithms introduced 
by Brown. Given the finite basis of such functions, it is extremely simple to determine any 
quantity that is defined by a power series expansion around the origin of the (w, w*) plane. Two 
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examples which we could evaluate with no ambiguity are the LL and NLL terms in the multi- 
Regge limit of the MHV amplitude. We could carry this exercise out through transcendental 
weight 10, and we presented the analytic formulae explicitly through six loops in Section HI The 
NMHV amplitudes also fit into the same mathematical framework, as we saw in Section [5j An 
integro-differential operator that generates the NMHV LLA terms from the MHV LLA ones jl3] 
has a very natural action on the SVHPLs, making it simple to generate NMHV LLA results to 
high order as well. A clear avenue for future investigation utilizing the SVHPLs is the NMHV 
six-point amplitude at next-to-leading-logarithm and beyond. 

A second thrust of this article was to understand the Fourier-Mellin transform from (w, w*) to 
(z/, n) variables. In practice, we constructed this map in the reverse direction: We built an ansatz 
out of various elements: harmonic sums and specific rational combinations of v and n. We then 
implemented the inverse Fourier-Mellin transform as a truncated sum, or power series around the 
origin of the (w,w*) plane, and matched to the basis of SVHPLs. We thereby identified specific 
combinations of the elements as building blocks from which to generate the full set of SVHPL 
Fourier-Mellin transforms. We have executed this procedure completely through weight six in the 
(z/, n) space, corresponding to weight seven in the (w, w*) space. In generalizing the procedure 
to yet higher weight, we expect the procedure to be much the same. Beginning with a linear 
combination of weight (p — 2) HPLs in a single variable x, perform a Mellin transformation to 
produce weight (p — 1) harmonic sums such as i/j, F 4 , F 6a , etc. For suitable combinations of these 
elements, the inverse Fourier-Mellin transform will generate weight p SVHPLs in the complex 
conjugate pair (w,w*). The step of determining which combinations of elements correspond to 
the SVHPLs was carried out empirically in this paper. It would be interesting to investigate 
further the mathematical properties of these building blocks. 

Using our understanding of the Fourier-Mellin transform, we could explicitly evaluate the 
NNLL MHV impact factor ^^(z/, n) which derives from a knowledge of the three-loop remainder 
function in the MRK limit (281 HO] • We then went on to four loops, using a computation of the 
four-loop symbol [22] in conjunction with additional constraints from the multi-Regge limit to 
determine the MRK symbol up to one free parameter a (which we suspect is zero). We matched 
this symbol to the symbols of the SVHPLs in order to determine the complete four-loop remainder 
function in MRK, up to a number of beyond-the-symbol constants. This data, in particular 
and then led to the NNLL BFKL eigenvalue E^n and N 3 LL impact factor <3?^| g (z/, n). These 
quantities also contain the various beyond-the-symbol constants. Clearly the higher-loop NNLL 
MRK terms can be determined just as we did at LL and NLL, using the master formula (12. 9 p 
and the SVHPL basis. However, it would also be worthwhile to understand what constraint 
can fix ao, and the host of beyond-the-symbol constants, since they will afflict all of these terms. 
This task may require backing away somewhat from the multi-Regge limit, or utilizing coproduct 
information in some way. 

We also remind the reader that we found that the NNLL BFKL eigenvalue E^X does not 
vanish as v — > 0, taking the limit after setting n = 0. This behavior is in contrast to what 
happens in the LL and NLL case. It also goes against the expectations in ref. [10], and thus calls 
for further study. 

Although the structure of QCD amplitudes in the multi-Regge limit is more complicated than 



51 



those of planar M = 4 super- Yang-Mills theory, one can still hope that the understanding of the 
Fourier-Mellin (u,n) space that we have developed here may prove useful in the QCD context. 

Finally, we remark that the SVHPLs are very likely to be applicable to another current 
problem in M = 4 super- Yang-Mills theory, namely the determination of correlation functions 
for four off-shell operators. Conformal invariance implies that these quantities depend on two 
separate cross ratios. The natural arguments of the poly logarithms that appear at low loop order, 
after a change of variables from the original cross ratios, are again a complex pair (w, w*) (or 
(z, z)). The same single- valued conditions apply here as well. For example, the one- loop off-shell 
box integral that enters the correlation function is proportional to L 2 (z, z)/(z — z). We expect 
that the SVHPL framework will allow great progress to be made in this arena, just as it has to 
the study of the multi-Regge limit. 
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A Single-valued harmonic polylogarithms 

A.l Expression of the L ± functions in terms of ordinary HPLs 

In this appendix we present the expressions for the Z 2 x Z 2 eigenfunctions L^(z) defined in 
eq. ( I3.19P as linear combinations of ordinary HPLs of the form H Wl (z) H W2 (z) up to weight 5. 
All expressions up to weight 6 are attached as ancillary files in computer- readable format. We 
give results only for the Lyndon words, as all other cases can be reduced to the latter. In the 
following, we use the condensed notation (I3.27P for the HPL arguments z and z to improve the 
readability of the formulas. 

A. 2 Lyndon words of weight 1 



L, 




o 




log |1 + w\ 2 + - log \w 



(A.l) 



(A.2) 



A. 3 Lyndon words of weight 2 



1 



1 



2 Hi q + 2 H\ q + 2 Hq H\ — 2 H Hi + 2 H 2 — 2 H 2 



(A.3) 



U 2 (z) - U 2 (z) + - log |z| 2 (log(l -z)- log(l - z)) , 
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A. 4 Lyndon words of weight 3 



-^2,1 



2 #0 H 0>0 + 2 # # 10 + 2 # # , + 2 # H 1>0 + 2 Hi # 0)0 +2 Hi # , (A. 4) 



+2 #0,0,0 + 2 #1,0,0 + 2 #0,0,0 + 2 #1,0,0 + 2 #3 + 2 # 3 



1 



Li 3 (z) + Li 3 (» - ~ log |2;| 2 Li 2 (z) + Li 2 (z) 



2 I „|2 



_ 4 lo S \ z 



log(l - z) + log(l - z) 



3 I „|2 



12 



log 1-2 



#0 #1,0 + #0 #1,0 + #1 #0,0 + #1 #0,0 + 2 #0 #0,0 + 2 #0 #1, 



+2 H # 0i o + 2 # #i,i + #1,0,0 + 2 #0,0,0 + 2 #2,0 + 2 # 2 ,i + 2 #1,1,0 
+#1,0,0 + 2 #0,0,0 + 2 #2,0 + 2 #2,1 + 2 #1,1,0 + 2 # #2 + 2 # #2 



+2#!# 2 + 2#! # 2 + #3 + #3 - 4Cs 

1 



-Li 3 (l - z) - Li 3 (l - 2) 



Li 3 (z +Li 3 (z 



+ ^og\z\ 2 



Li 2 (z) + Li 2 (z) 



--log|l-z| 2 



Li 2 (z) + Li 2 (z) 



^log|z| 2 log|l-z| 2 + ^log 3 |z| 2 



■TlOg- 

4 z 



log 2 (l - z) - log 2 (l - z)] + C2 log |1 - z\ 2 + C; 



3 • 



A. 5 Lyndon words of weight 4 



(A.5) 



Lti 



#0 #2,0 + #0 #1,0,0 — #0 #2,0 — #0 #1,0,0 — #1 #0,0,0 + #1 #0,0,0 



(A.6) 



+#0,0 #2 + #0,0 #1,0 — #0,0 #2 — #0,0 #1,0 + 2 #0 #1,1,0 — 2 #0 #1,1,0 
+2 #0,0 #1,1 — 2 #0,0 #1,1 + #3,0 ~~ #2,0,0 ~~ #1,0,0,0 + 2 #3,1 — 2 #1,1,0,0 
+#2,0,0 + #1,0,0,0 — #3,0 — 2 #3,1 + 2 #1,1,0,0 — #0 #3 + #0 #3 — 2 #1 #3 

+2 H 1 # 3 + 4 #1 Cs + #4 - 4 #! Cs " #4 



2 #0 #1,0,0 — 2 #0 #1,0,0 — 2 #1 #0,0,0 + 2 #1 #0,0,0 + 2 #0,0 #1,0 



(A.7) 



-2 #0,0 #1,0 — 2 #1,0,0,0 + 2 #1,0,0,0 + 2 # 4 — 2 # 4 
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^2,1,1 



#o #1,0,0 + #o Hi 2 + #o #i,i,o — #o #1,0,0 — #0 #1,2 — #0 #1,1,0 



(A. 



-#! #0,0,0 — #1 #2,0 + #1 #0,0,0 + #1 #2,0 + #0,0 #1,0 + #0,0 #1,1 
"#0,0 #1,0 — #0,0 #1,1 + #2 #1,0 — #2 #1,0 + 2 #0 #1,1,1 — 2 #0 #1,1,1 

-2 #1 #2,1 + 2 #1 #2,1 + 2 #2 #1,1 — 2 #2 #1,1 + #3,1 + #2,2 

#1 ,0,0,0 — #1,2,0 — #1,1,0,0 + 2 #2,1,1 ~ 2 #1,1,1,0 + #1,0,0,0 + #1,2,0 + #1,1,0,0 
'#3,1 — #2,2 — 2 #2,1,1 + 2 #1,1,1,0 — #1 #3 + #1 #3 + 2 #1 (3 + #4 

-2 #1 C3 — #4 



A. 6 Lyndon words of weight 5 



Li 



T + 

-^3,1, 



2 #0 #0,0,0,0 + 2 #0 #1,0,0,0 + 2 #0 #0,0,0,0 + 2 #0 #1 



0,0,0 



2 #1 #o,o,o,( 



(A.9) 



+2#!# ,0,0,0 + 2 #0,0 #0,0,0 + 2 #0,0 #1,0,0 + 2 #0,0 #0,0,0 + 2 #0,0 #1,0,0 

+2 #1,0 #0,0,0 + 2 #1,0 #0,0,0 + 2 #0,0,0,0,0 + 2 #1,0,0,0,0 + 2 #0,0,0,0,0 + 2 #1,0,0,0,0 



+2 # 5 + 2 # £ 



#5 + #5 + #4,0 + #4,0 + #4,1 + #4,1 + #3,2 + #3,2 + #3,1, 



# 



3,1,0 



(A.10) 



-# 



2,0,0,0 



#2,0,0,0 + #2,1,0,0 + #2,1,0,0 + #1,0,0,0,0 + #1,0, 



0,0,0 



#1 



,2,0,0 



+#1,2,0,0 + #1,1,0,0,0 + #1,1,0,0,0 + 2 #0,0,0,0,0 + 2 #0,0,0,0,0 + 2 #3,0,0 + 2 #3,0,0 
+2 #3,1,1 + 2 #3,1,1 + 2 #1,1,1,0,0 + 2 #1,1,1,0,0 + 4 C5 + #0 #4 + #0 #3,1 + #0 #2,0,0 
+#0 #2,1,0 + #0 #1,0,0,0 + #0 #1,2,0 + #0 #1,1,0,0 + #0 #4 + #0 #3,1 + #0 #2,0,0 
+#0 #2,1,0 + #0 #1,0,0,0 + #0 #1,2,0 + #0 #1,1,0,0 + #1 #0,0,0,0 + #1 #4 + #1 #3,0 
+#1 #0,0,0,0 + #1 #4 + #1 #3,0 + #0,0 #2,0 + #0,0 #2,1 + #0,0 #1,0,0 + #0,0 #1,2 
+#0,0 #1,1,0 + #0,0 #2,0 + #0,0 #2,1 + #0,0 #1,0,0 + #0,0 #1,2 + #0,0 #1,1,0 
+#2 #0,0,0 + #2 #3 + #2 #0,0,0 + #2 #3 + #1,0 #0,0,0 + #1,0 #3 + #1,0 #0,0,0 
+#1,0 #3 + #1,1 #0,0,0 + #1,1 #0,0,0 + 2 #0 #0,0,0,0 + 2 #0 #3,0 + 2 #0 #1,1,1,0 
+2 # #0,0,0,0 + 2 #0 #3,0 + 2 # #1,1,1,0 + 2 #1 #3,1 + 2 #! # 3;1 + 2 # 0i0 #0,0,0 
+2 #0,0 #3 + 2 #0,0 #1,1,1 + 2 #0,0 #0,0,0 + 2 #0,0 #3 + 2 #0,0 #1,1,1 — 2 # 2 C3 
-2 # 2 Cs - 2 #1,0 Cs - 2 #1,0 C 3 + 2 #1,1 #3 + 2 #!,i # 3 - 4 # 1A C 3 - 4 # w C 3 
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—2 H H i ( 3 — 2 H Hi ( 3 



#5 + #5 + # 4 ,i + # 4 ,i + #2,3 + #2,3 + #1,0,0,0,0 + #1,0,0,0,0 + #1 



3,0 



(A.11) 



+#1,3,0 + #1,1,0,0,0 + #1,1,0,0,0 + 2 #0,0,0,0,0 + 2 #0,0,0,0,0 + 2 #4,0 + 2 i? 4i o 

+2 #2,0,0,0 + 2 #2,0,0,0 + 2 #2,2,0 + 2 #2,2,0 + 2 #2,2,1 + 2 #2,2,1 + 2 #1,1,2,0 

+2 #1,1,2,0 — 6 (5 + #0 #1,0,0,0 + #0 #1,3 + #0 #1,1,0,0 + #0 #1,0,0,0 + #0 #1,3 

+#0 #1,1,0,0 + #1 #0,0,0,0 + #1 #4 + #1 #2,0,0 + #1 #0,0,0,0 + #1 #4 + #1 #2,0,0 

+#0,0 #1,0,0 + #0,0 #1,1,0 + #0,0 #1,0,0 + #0,0 #1,1,0 + #2 #1,0,0 + #2 #1,0,0 

+#1,0 #0,0,0 + #1,0 #2,0 + #1,0 #0,0,0 + #1,0 #2,0 + #1,1 #0,0,0 + #1,1 #0,0,0 

+2 # #0,0,0,0 + 2 # #4 + 2 # #2,0,0 + 2 H # 2 , 2 + 2 H #1,1,2 + 2 if #0,0,0,0 

+2 #0 #4 + 2 #0 #2,0,0 + 2 #0 #2,2 + 2 #0 #1,1,2 + 2 #1 #2,2 + 2 #1 #2,2 

+2 #0,0 #0,0,0 + 2 #0,0 #2,0 + 2 #0,0 #0,0,0 + 2 #0,0 #2,0 + 2 #2 #0,0,0 + 2 #2 #2,0 

+2 # 2 #1,1,0 + 2 #2 C3 + 2 #2 #0,0,0 + 2H 2 #2,0 + 2 # 2 #1,1,0 + 2 #2 C3 

+2 #1,1 #2,0 + 2 #1,1 #2,0 — 4 # 00 ( 3 — 4 # , ( 3 + 4 # lj0 C3 + 4 #1,0 C3 



-8 # 1 C3 + 8 #11 C3 — 4 H H ( 3 + 4 # #1 C3 + 4 #0 #1 C3 



#0 #2,0,0 + #0 #1,0,0,0 + #0 # 



-"2,0,0 



#0 #1,0,0,0 + #1 #0. 



0,0,0 



#1 #0,0,0,0 (A. 12) 



+#0,0 #2,0 + #0,0 #1,0,0 + #0,0 #2,0 + #0,0 #1,0,0 + #2 #0,0,0 + #2 #0,0,0 
+#1,0 #0,0,0 + #1,0 #0,0,0 + 2 Hq #3,0,0,0 + 2 Ho #1,1,0,0 + 2 #0 #0,0,0,0 
+2 #0 #1,1,0,0 + 2 #0,0 #0,0,0 + 2 #0,0 #1,1,0 + 2 #0,0 #0,0,0 + 2 #0,0 #1,1,0 
+2 #1,1 #0,0,0 + 2 -f/"i,i #0,0,0 — 4 #0,0 C3 — 4 #i,o C3 + #4,0 + #2,0,0,0 
+#1,0,0,0,0 + 2 #0,0,0,0,0 + 2 #4,1 + 2 #1,1,0,0,0 — 4 #0,0 C3 — 4 #1,0 C3 + #4,0 



+# 



2,0,0,0 



#1 



0,0,0,0 



2#, 



0,0,0,0,0 



2 #4,1 + 2 #1,1,0,0,0 — 4 # # C3 



-4#o#iC3-4#o#lC3 + #0#4 + #0#4 + 2#i#4 + 2#i# 4 
-# 5 + #5-4(5 



#0 #1 



0,0,0 



#0 #1, 



0,0,0 



#!# 



0,0,0,0 



#1#, 



0,0,0,0 



#0,0 #1,0,1 



(A.13) 



+#0,0 #1,0,0 + #1,0 #0,0,0 + #1,0 #0,0,0 + 2 #0 #0,0,0,0 + 2 #0 #3,0 + 2 H #1,2,0 
+2 #0 #0,0,0,0 + 2 #0 #3,0 + 2 Ho #1,2,0 + 2 #1 #3,0 + 2 Hi # 3 , + 2 # ,o #0,0,0 
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+2 # , #3 + 2 # 0)0 #1,2 + 2 #o,o #o,o,o + 2 #0,0 #3 + 2 # 0>0 #1,2 + 2 #i, # 3 
+2 #1,0 #3 + 8 #0,0 C3 + 8 #1,0 C3 + #1,0,0,0,0 + 2 #0,0,0,0,0 + 2 #3,0,0 + 2 #3,2 
+2 #1,2,0,0 + 8 #0,0 C3 + 8 #1,0 C3 + #1,0,0,0,0 + 2 #0,0,0,0,0 + 2 #3,0,0 + 2 #3,2 



+ 2#x 



,2,0,0 



8# #oC3 + 8#o#lC3 + 8#o#lC3 + #5 + #5 + 16C E 



#5 + #5 + #4,1 + #4,1 + #3,2 + #3,2 + #3,1,1 + #3,1,1 + #2,3 + #2,: 



(A.14) 



+#2,2,1 + #2,2,1 + #2,1,2 + #2,1,2 + #1,0,0,0,0 + #1,0,0,0,0 + #1,3,0 + #1,3,0 
+#1,2,0,0 + #1,2,0,0 + #1,2,1,0 + #1,2,1,0 + #1,1,0,0,0 + #1,1,0,0,0 + #1,1,2,0 + #1,1,2,0 

+#1 ,1,1,0,0 + #1,1,1,0,0 + 2 #0,0,0,0,0 + 2 #0,0,0,0,0 + 2 #4,0 + 2 # 4i0 + 2 #3,0,0 

+2 #3,0,0 + 2 #3,1,0 + 2 #3,1,0 + 2 #2,0,0,0 + 2 #2,0,0,0 + 2 #2,2,0 + 2 #2,2,0 

+2 #2,1,0,0 + 2 #2,1,0,0 + 2 #2,1,1,0 + 2 #2,1,1,0 + 2 #2,1,1,1 + 2 #2,1,1,1 + 2 #1,1,1,1,0 

+2 #1,1,1,1,0 — 4 £5 + #0 #1,0,0,0 + #0 #1,3 + #0 #1,2,0 + #0 #1,2,1 + #0 #1,1,0,0 

+#0 #1,1,2 + #0 #1,1,1,0 + #0 #1,0,0,0 + #0 #1,3 + #0 #1,2,0 + #0 #1,2,1 

+#0 #1,1,0,0 + #0 #1,1,2 + #0 #1,1,1,0 + #1 #0,0,0,0 + #1 #4 + #1 #3,0 + #1 #3,1 

+#1 #2,0,0 + #1 #2,2 + #1 #2,1,0 + #1 #0,0,0,0 + #1 #4 + #1 #3,0 + #1 #3,1 

+#1 #2,0,0 + #1 #2,2 + #1 #2,1,0 + #0,0 #1,0,0 + #0,0 #1,2 + #0,0 #1,1,0 

+#0,0 #1,1,1 + #0,0 #1,0,0 + #0,0 #1,2 + #0,0 #1,1,0 + #0,0 #1,1,1 + #2 #1,0,0 

+#2 #1,2 + #2 #1,1,0 + #2 #1,0,0 + #2 #1,2 + #2 #1,1,0 + #1,0 #0,0,0 + #1,0 #3 

+#1,0 #2,0 + #1,0 #2,1 + #1,0 #0,0,0 + #1,0 #3 + #1,0 #2,0 + #1,0 #2,1 + #1,1 #0,0,0 

+#1,1 #3 + #1,1 #2,0 + #1,1 #0,0,0 + #1,1 #3 + #1,1 #2,0 + 2 #0 #0,0,0,0 + 2 #0 #4 

+2 # #3,0 + 2 # #3,1 + 2 # #2,0,0 + 2 #0 #2,2 + 2 #0 #2,1,0 + 2 #0 #2,1,1 

+2 # #1,1,1,1 + 2 # #0,0,0,0 + 2 #0 #4 + 2 # #3,0 + 2 # #3,1 + 2 # #2,0,0 

+2 # #2,2 + 2 #0 #2,1,0 + 2 #0 #2,1,1 + 2 #0 #1,1,1,1 + 2 #1 #2,1,1 + 2 #1 #2,1,1 

+2 #0,0 #0,0,0 + 2 #0,0 #3 + 2 #0,0 #2,0 + 2 #0,0 #2,1 + 2 #0,0 #0,0,0 + 2 #0,0 #3 

+2 #0,0 #2,0 + 2 #0,0 #2,1 + 2 #2 #0,0,0 + 2 #2 #3 + 2 #2 #2,0 + 2 #2 #2,1 

+2 #2 #1,1,1 + 2 #2 #0,0,0 + 2 #2 #3 + 2 #2 #2,0 + 2 #2 #2,1 + 2 #2 #1,1,1 

+2 #1,1 #2,1 — 2 #1,1 (3 + 2 #1,1 #2,1 — 2 #1,1 £3 



56 



A. 7 Expression of Brown's SVHPLs in terms of the L ± functions 

In this appendix we present the expression of Brown's SVHPLs corresponding to Lyndon words 
in terms of the Z2 x Z2 eigenfunctions L^(z). 

£0 — L Q , 

L -L+ V 
^1 — % — 2 ' 

L 2 = L 2 , 

£3 = Lt - — [Lq] 3 , 

£2,1 = L\ [Lq ] 2 + - £3 + L 21 + (3 , 

La — £4 , 
£3,1 = — ^ £2 [£q ] 2 + L 4 + , 



£2,1,1 — — ^ £2 L Lf + - + £2,1,1 
1 



L5 ~ L ^"2i0 [Lo]5 ' 
£4,1 = [£o ] 4 " J £ 3 + [^0 ] 2 + 5 [L f Cs + § ^ + ^1 + C 5 , 

£3,2 = ~ Lf [L -] 4 + ± Li [L-f -I Li- [L- ] 2 Cs + ^2 - 4Cb , (A ' 15) 

£3,1,1 = Tjj [£o ] 3 [£i"] 2 _ ^ £2,1 [Lq] 2 + — [£o ] 5 - ^Lq Li Li + -L Q Li ( 3 
+ £4,1 + £3^1,1 > 

£2,2,1 = -— [Lq] 3 [Li] 2 + - L^x [Lq] 2 - [Lq] 5 + - Lq Li Li - -L Li ( 3 

- 7 L~ - 1 L~ + L+ 

2 -^4,1 2 3 ' 2 2,2,1 5 

£2,1,1,1 = — [Lq] 2 [Li] 3 - — £1" [£q ] 4 + ^ £3 [£o ] 2 - g [£o ] 2 Cs - y £0 £2,1 £| 

— 7 £5 + 2 ^1,1 + 2 Cs + £2,1,1,1 j 
Lq = L 6 , 

£5,1 = - J £4" [£o ] 2 + ^ £ 2 " [£o ] 4 + 2 £ 6 ~ + L+! , 



£4,2 -7/4 [£q ] 2 ~ Y2" £2 [£q ] 4 2~ + ^ 2 ' 
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1 3 1 1 1 3 

£4,1,1 = —L 2 L 1 [L } — - L 31 [L \ — - L 4 L Q L 1 + - L 2 L ( 3 + - L 51 + £ 411 , 

3 3 1 2 3 1 7 

£3,2,1 = -~^L 2 L 1 [L ] + - £ 31 [L ] + - L 4 L Q L 1 — - L 2 L Q ( 3 - - L 51 + L 321 , 

1 3 

£3,1,2 = — ^ ^2 A) £3 — 2 ^ 2 ^° ^ 3 + ^ 3 ' 1 ' 2 + ^ ^S' 1 + "^ 4 > 2 ' 

£3,1,1,1 =J^ L 2 l L o ] 2 + - L 4 [L ] 2 - — L 2 [L Q ] 4 - - L 2)1)1 [£ ] 2 - -L L\ L\ x 

— L 6 + £4,1,1 + ^3,1,1,1 > 

1 2 2 ^ 2 1 4 ^ 2-^1 

£2,2,1,1 = — ^ £2 [£q] [£l] _ ^£4[£o] +Y2^ 2 ^°] +4 ^2,1,1 [£ol +^^6 

+ 7 £2 £l" £3 — o ^ 2 ^l" ^ 3 ~*~ 4 ^ ^1 ~~ 2 ^ 3 ' 1,2 ~~ ^ ^M.l ~~ -^3,2,1 + £^2,1,1 J 

£2,1,1,1,1 = — £2 £1" [£o ] 3 + yjj £ji [£o ] 2 + 4^ £2 £0 [£f] 3 + g £4 £0 £^ _ 4 £2 £0 C3 

— 4 £0 £2,1,1 £1" ~~ 4 £5^1 + 2 ^1,1,1 + £2,1,1,1,1 • 

(A.16) 

B Analytic continuation of harmonic sums 

In this section we review the analytic continuation of multiple harmonic sums and the structural 
relations between them, as presented by Bliimlein [50]. Multiple harmonic sums are defined by, 

N ki fcn-l / \l, / \k 

a (m STX^ ^W^iT sgn(a w ) " 

fc 1= l fc 2 =l fe„=l K l K n 

where the at are positive or negative integers, and N is a positive integer. For the cases in which 
we are interested, they are similar to the Euler-Zagier sums (I3.10p . except that the summation 
range differs slightly. They are related to Mellin transforms of real functions or distributions 



S ai ,...,a n (N) 



C dxx N f ai _ an =M[f ai _ an (x)](N). (B.2) 
Jo 



Typically f(x) are HPLs weighted by factors of 1/(1 ± x). To avoid singularities at x — 1, it is 
often useful to consider the +-distribution, 

M[(f(x)) + ](N) = f dx {x N - 1) f(x) . (B.3) 



The weight \w\ of the harmonic sum is given by \w\ = ^2=1 \ a k\- The number of harmonic sums 
of weight w is equal to 2 • 3' u, ' _1 , but not all of them are independent. For example, they obey 
shuffle relations [72J. It is natural to ask whether these are the only relations they satisfy. In fact, 
it is known that in the special case N — > oo, in which the sums reduce to multiple zeta values, 



58 



many new relations emerge [731 EH ED ES]- In ref. [50], an analytic continuation of the harmonic 
sums was considered. It is defined by the integral respresentation, eq. (IB.2p . where N is allowed 
to take complex values. This allows for two new operations — differentiation and evaluation at 
fractional arguments — which generate new structural relations among the harmonic sums. 

In the present work, harmonic sums with negative indices do not appear, so we will assume 
that cik > 0. This assumption provides a considerable simplification. The derivative relations 
allow for the extraction of logarithmic factors, 



M[k>g'(x) f(x)](N) 



dN l 



M[f(x)](N) 



(B.4) 



which explains why the derivatives of the building blocks in Section |6] generate SVHPLs. In 
ref. [50] , all available relations are imposed, and the following are the irreducible functions through 
weight five: 



weight 1 



weight 3 



weight 4 



S 1 (N)=^j(N + 1)+ 1e = M 



F 4 (N) = M 



F 6a (N) = M 
F 7 (N) = M 



x — 1 



(TV) 



Li 2 ( 


z)\ 


1 - 


X ) \ 


Li 3 ( 


x)\ 




X J . 


£l,2 





(N) 



(B.5) 



:b.6) 



x 



(N) 
(N) 



(B.7) 



weight 5 



F g (N) = M 
F n (N) 
F 13 (N) 
F 17 (N) 



M 



M 



M 



Li 4 (a;) 



x — 


1 J 


£2,2 




x - 




ul{ 




x — 




^1,3 


(x)\ 


x - 





(N) 
(N) 
(N) 
(N) 



(B.8) 



There are no irreducible basis functions of weight two. These functions are meromorphic with 
poles at the negative integers. To use these functions in the integral transform (14. 4p . we need 
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the expansions near the poles. Actually, we only need the expansions around zero, since the 
expansions around any integer can be obtained from them using the recursion relations of ref. 



^ n \l + z 
F 4 (z 
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4z 2z 3 
Fi 3 (z — 1) H ^ 5-^1(2; 



2S 2 Az) , 2 



2 



F 17 (z-1) 



Ci 



[5 1 3 (^) + 35x^)52^) +25 3 (z)] 



2 6z 2 

The expansions around zero can be obtained from the integral reprsesentations. We find that, 
for 5 — > 0, the expansions can all be expressed simply in terms of multiple zeta values, 
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n=l 
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These single- variable functions can be assembled to form two-variable functions of v and n, such 
that their inverse Mellin-Fourier transforms produce sums of SVHPLs. This construction is not 
unique, because other building blocks could be added. We choose to define the two-variable 
functions as, 
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B.l The basis in (z/, n) space in terms of single-valued HPLs 

In this appendix we present the analytic expressions for the basis of Z 2 x Z 2 eigenfunctions in 
(u, n) space in terms of single-valued HPLs in (w,w*) space up to weight five. The Z 2 x Z 2 
acts on (w,w*) space via inversion, while it acts on {v,n) space via via [n -H- — n, v —v\ and 
[v -H- — 1/]. The eigenvalue under Z 2 x Z 2 in (w, w*) space will be referred to as parity. 

Basis of weight 1 with parity (+,+): 



X[l] = 2L+. 



(B.13) 



Basis of weight 1 with parity (+, — ): 



X [<5 ,n] 



L -. 



(B.14) 
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Basis of weight 2 with parity (+,+): 

X[E v , n ] = [L+] 2 -i[L "] 2 , (B.15) 
X[6 , n /(iv)} = l -[L f. (B.16) 

Basis of weight 2 with parity (+,—): 

X[V] =-L L+. (B.17) 

Basis of weight 2 with parity ( — ,—): 

X[N] =4^2 . (B.18) 

Basis of weight 3 with parity (+,+): 

IKJ = l[Lt?-Lt, (B.19) 
X[N 2 ] = 12 Lt - 2 Lf[L ]\ (B.20) 

X[V 2 } = ~ Lt [L -] 2 - Lf . (B.21) 

Basis of weight 3 with parity (+,—): 

X[VE v>n ] = i[Lo] 3 -2L^, (B.22) 
X[D v E u , n ] = -^[L ] 3 -L -[L+] 2 + 4L 2 - 1 , (B.23) 
X [5 , n /M 2 ] = i[L ] 3 . (B.24) 

Basis of weight 3 with parity (—,+): 

X[NV] = -LJ Lq . (B.25) 

Basis of weight 3 with parity (—,—): 

Z[iV = 2^2 L+ . (B.26) 
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Basis of weight 4 with parity (+,+): 
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Basis of weight 5 with parity (+, + 
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Basis of weight 5 with parity 
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Basis of weight 5 with parity ( 
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